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Abstract 

Consider a system of N bosons on the three dimensional unit torus interacting via a pair 
potential N'^V{N{xi—Xj)), where x = [xi, . . . , xn) denotes the positions of the particles. Suppose 
that the initial data il^Nfi satisfies the condition 

where Hjv is the Hamiltonian of the Bose system. This condition is satisfied if ipNfi — Wm4>n.o 
where Wn is an approximate ground state to Hf^ and (pN.o is regular. Let ^N,t denote the 
solution to the Schrodinger equation with Hamiltonian Hj^. Gross and Pitaevskii proposed to 
model the dynamics of such system by a nonlinear Schrodinger equation, the Gross-Pitaevskii 
(GP) equation. The GP hierarchy is an infinite BBGKY hierarchy of equations so that if Ut 
solves the GP equation, then the family of fc-particle density matrices {®kUt:k > 1} solves the 
GP hierarchy. We prove that as — > cx) the limit points of the /c-particle density matrices of 
ipN,t are solutions of the GP hierarchy. The uniqueness of the solutions to this hierarchy remains 
an open question. Our analysis requires that the N boson dynamics is described by a modified 
Hamiltonian which cuts off the pair interactions whenever at least three particles come into a 
region with diameter much smaller than the typical inter-particle distance. Our proof can be 
extended to a modified Hamiltonian which only forbids at least n particles from coming close 
together, for any fixed n. 



1 Introduction 

A very simple and useful way to understand Bose systems is to treat all bosons as independent 
particles. In particular, to prove the Bose-Einstein condensation is a simple exercise in the case of 
non-interacting bosons J2]- The true many-body problem with a pair interaction is a much harder 
problem. Gross [HI El and Pitaevskii |^ proposed to model the many-body effect by a nonlinear 
on-site self interaction of a complex order parameter (the "condensate wave function" ) . The strength 
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of the nonlinear interaction in this model is given by the scattering length ag of the pair potential. 
The Gross-Pitaevskii (GP) equation is given by 



■« A , I i2 S£{u,u) 

idtut = -Aut + a\ut\ ut ^ 



5u 



£{u,u) = j dx \Vu\'^ + ^\u\^ , (1.1 



where £ is the Gross-Pitaevskii energy functional and a = Syrao. The Gross-Pitaevskii equation was 
considered as a mean-field model; some corrections were obtained in |14[ ll5j and more recently, e.g., 
in [1H|. 

The many-body effects can be analyzed in much more details by the Bogoliubov transforma- 
tion. The Bogoliubov's theory, it should be emphasized, postulates that the ratio between the 
non-condensate and the condensate is small. Apart from this assumption, the coupling constant a 
derived from the Bogoliubov's theory is the semiclassical approximation to the scattering length. 
One can perform some higher order diagrammatic re-summation to recover the scattering length. 
In the case of hard core potential, the traditional theory relies on the non-rigorous pseudo-potential 
model 13 . 

The first rigorous result concerning the many-body effects of the Bose gas was Dyson's estimate 
of the ground state energy. Dyson proved the correct leading upper bound to the energy and a 
lower bound off by a factor around 10. The matching lower bound was obtained by Lieb-Yngvason 
[241 125j forty years later. The last result has inspired many subsequent works, including a proof 
[13 EH Eg that the Gross-Pitaevskii energy functional correctly describes the ground state in a 
scaling limit to be specified later. 

The experiments on the Bose-Einstein condensation were conducted by observing the dynamics of 
the condensate when the confining traps are switched off. It is remarkable that the Gross-Pitaevskii 
equation, despite being a mean-field equation, has provided a very good description for the dynam- 
ics of the condensate. The validity of the Gross-Pitaevskii equation asserts that the fundamental 
assumption of the Gross-Pitaevskii theory (i.e., that the many-body effects can be modelled by a 
nonlinear on-site self interaction of the order parameter) applies not only to the ground states, but 
to certain excited states and their subsequent time evolution. This remarkable and fundamental 
property of the Bose gas has mostly been taken for granted and has not been treated rigorously in 
the literature. In order to explain the issues involved, we now introduce some notations and set up 
the scaling for the Gross-Pitaevskii theory. 

Let y > be a fixed nonnegative, spherically symmetric, smooth potential with compact support. 
The zero energy scattering solution / satisfies the equation 

A + ^V{x)]f{x) = 0. (1.2) 

If we fix the normalization lim|a,|^oo fix) = 1 and write /(r) = q(r)/r (where r = \x\), the scattering 
length Oo of V is defined by 

ao = lim r — q{r) (1-3) 

r— >oo 

and thus for x large, 

/(x)~l-ao/|x|. (1.4) 
From the zero energy equation, we also have the identity 

dxV{x)f{x) = Svrao . (1.5) 

Let A be the three-dimensional torus of unit side-length. Denote the position of bosons in A 
by X = {xi,X2, ■ ■ ■ ,xn), Xj E A. Lieb, Seiringer and Yngvason ^\ pointed out that, in order to 
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obtain the GP functional, the length scale of the pair potential should be of order Notice that 

the density of the system is N and the typical inter-particle distance is N^^^'^, which is much bigger 
than the length scale of the potential. The system is really a dilute gas on the scale of the range of 
the interaction, but it is scaled in such a way that the size of the total system is of order one. 
The Hamiltonian of the Bose system on the torus A is given by 

N N 

= - J^A, + J^iVV(iV(x, -Xfc)) . (1.6) 
j=i j<k 

By scaling, the scattering length of the potential N'^V{Nx) is a := uq/N. We shall from now on use 
the notation 

Va{x) ■.= N^V{Nx) , a = ao/N. (1.7) 
Notice that the pair potential in ()1.6|) is an approximation to the mean field Dirac delta interaction: 

b ^ r 

-^^5{xj-Xk), bo= dxV{x). 

j<k 

The Schrodinger equation is given by 

idtilJt = Hipt, or idtjN = [H,-fN] , (1-8) 

where 'j^ is the A^-particle density matrix. For a pure state, 7^ := is the orthogonal projection 

onto ijj. 

Introduce the shorthand notation 

X := {xi,X2,. . . ,xn), Xfc := (xi, . . . ,rEfc), XN-k ■= {xk+i, ■ ■ ■ ,xn) 
and similarly for the primed variables, x'^ := (x'l, . . . , x^). The /c— particle density matrix is given by 



(fc) 



{^k;^'k) ■= j 'i^N-kiN {^k,^N-k;^'k,^N-k) ■ (1-^ 



Our normalization implies that Tr = 1 for all k = 1, . . . , N. 
The density matrix 7^^^ satisfies the following equation 

idt-f'-ji\{xi; x[) = (-A^.j + A^>^)-f^j^\{xi; x[) 

f (2) (1-10) 

+ (iV - 1) / dX2{Va{xi - X2) - Va{x[ - X2))TN^^{xl, X2; x[, X2) , 

and similar equations hold for 7^^ for A; > 1. To close this equation, one needs to assume some 
relation between 7^'*^ and 7^'*^. The simplest assumption would be the factorization property, i.e., 

j^^l{xi,X2;x[,x'2) = 7£(a;i;a;i)7^^^,(x2;x2) . (1.11) 

(k) (k) 

This does not hold for finite N, but it may hold for a limit point 7^ of 7]^ ^ as ^ cxd, i.e., 

7p^(xi,X2;x'i,X2) = 7^ ^(xi;x'i)7f ^(x2;x2) . (1.12) 
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Under this assumption, 7^^^^ satisfies the hmiting equation 

idat^\xi;x[) = (-A^i + A^0j^^\xi;x[) + {Qt{xi) - Qt{x'i))j^^\xi; x[) (1.13) 

where 

Qt{x) = lim N [ dyVaix-y)ptiy), pt{x) = ^i^\x; x) . (1.14) 

If, instead of V^(x), we use the unsealed mean field potential V{x)/N, then Qt is the convolution 
of V with the density pt{x). The equation H1.13() becomes the Hartree equation. For Bose systems 
with mean field interaction and product initial wave function, the factorization assumption H1.12() 
can be proved for a general class of potentials. See the work of Hepp and Spohn j21j for bounded 
potentials and |S] for potentials with Coulomb type singularity. For certain quasi-free initial data, 
Ginibre and Velo ^ can handle all integrable singularities. We note that in one dimension the 
convergence to the GP hierarchy 1)1.24(1 for the delta potential was established by Adami, Bardos, 
Golse and Teta in pi. 

In our setting, Va{x) lives on scale 1/A^ and is much more singular than all the cases considered 
previously. If pt{x) is continuous, then Qt is given by 

Qt{x) = bopt{x). 

Thus (|1.13() gives the GP equation with the incorrect coupling constant cr = 60 instead of o" = Svrao. 
It is known that bo /Sir is the first Born approximation to the scattering length oq and the following 
inequality holds: 

ao<^ = ^ [ lv{x)dx. (1.15) 
oTT 47r Jr3 2 

To go beyond the Born approximation, we need to understand the short scale correlations of the 
ground state which we now review. 

The ground state of a dilute Bose system with interaction potential Va is believed to be very 
close to the form 

W{^):=l[f{Nix,-x,)) (1.16) 

i<j 

where / is the zero-energy solution (|1.2() . We remark that Dyson [S] took a different function which 
was not symmetric, but the short distance behavior was the same as in W. Since in the experiments 
the initial states were prepared with a trapping potential, living on a scale of order one, the ground 
state for such a trapped gas is of the form tpix) = W{x)(j)(x) ^1120] where (/>(x) is close to a product 
function. Thus we shall consider initial data of the form W{x)(p(x). 

We assume for the moment that the ansatz, ipt{^) = W{x)(pt{x.) with (pt a product function, 
holds for all time. The reduced density matrices for il^ti^) satisfy 

ji^\xi,X2-,x[,x'^)^f{N{x,-X2))f{N{x[-x'^)hl'\xi-,x[hl'\x2;x'^). (1.17) 
Together with (|1.5|) and the assumption that pt is smooth, we have 

lim N / dx2Va{xi - X2h''^\{xi,X2;x[,X2) = 8'n:ao-ft^\xi;x[)pt{xi) . (1.18) 

N^oo J ' 

We have used that lim|2,|^oo /(x) = 1 and the last equation is valid for |xi — x[\ » 1/A^. This gives 
the GP equation with the correct dependence on the scattering length. 

Notice that the relations (jl.l7|) and (jl.l8|) are very subtle. The correlation in 7'-^) occurs at the 
scale 1/A^. Testing the relation (|1.17() in a weak limit, all correlations at the scale 1/N disappear and 
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the product relation ()1.12() will hold. However, this short distance correlation shows up in the GP 
equation due to the singular potential NVa{xi — X2)- Therefore, a rigorous justification of the GP 
equation requires a proof that the relation (|1.17jl holds with such a precision that H1.18|) is valid — a 
formidable task. 

We would like to divide this task into two parts. The first question is whether the short scale 

(2) 

structure of is given by f{N{xi — X2))f{N{x'i — x'2)), i.e., whether, for large A^, 

7^^^j(a;i,X2;xi,4) ~ f{N{xi -X2))f{N{x^ - X2))7p''(xi, X2; x'^, Xg) (1.19) 

where 7^ is a weak limit of t i^^^ short distance correlations given by f{N{xi — X2))f{N{x'i — x'2)) 
vanish when the weak limit is taken) living on a scale of order one so that 

dx2V;(xi - X2)7]vj(xi,X2;xi,X2) = 87rao7t \xi,xi\x'i,xi) . (1.20) 

The second question is whether 7^^ factorizes, i.e., whether 

lf\xi,X2]Xi,X2) = 7^ ^(xi;x'i)7f ^(X2;x2). (1.21) 

If we replace the Hamiltonian (|1.6() by a modified Hamiltonian in which we remove the pair 
interaction when three or more particles come close together in a very short distance, then we can 
prove a certain version of (|1.19|) and (|1.2fl|) . Let i denote a distance much smaller than the typical 
inter-particle distance, say, I = N^^^'^^^ for some 6 > 0. The modified Hamiltonian is approximately 
of the form 

~ - ^ A,- + ^ [n^V{N{xj - Xk)) n - Xk\ > £)] . (1.22) 

j=l j<ck i^j,k 

The precise definition will be given in (|2.9|) . This cutoff modification changes the Hamiltonian for 
events that are rare with respect to the expected typical particle distribution, therefore it should have 
little effect on the dynamics. Unfortunately we cannot control this effect rigorously. (In principle, the 
original unmodified dynamics may introduce local clustering of particles, despite that it is unfavorable 
for the local energy.) In the computation of the ground state energy by Lieb-Yngvason 24 , no such 
modification was needed since the positivity of the contribution from these rare events could be 
exploited. Our method is based on the conservation of along the dynamics and an inequality of 
the form: 

Y, j W^\V.iVj(l,t\^ <C I \{H + N)i;t\^ = C J \{H + N)4^o\^ < C N\ (1.23) 

i<j 

with V't = W(j)f The idea of using the conservation of higher power of the Hamiltonian was introduced 
in [nj. Clearly, the computation of involves derivatives of the pair potential which have no definite 
sign. If we use the original Hamiltonian H instead of H, these terms cannot be controlled by the 
kinetic energy operator in the rare situation when many particles come close together. The modified 
Hamiltonian (|1.22|) removes this technical obstacle. 

Using (|1.23|) we will prove that weak limit points 7^^^"* of the /c-particle density 7^^ (whose 
time evolution is generated by the modified Hamiltonian H) satisfy the following infinite BBGKY 
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hierarchy, which will be called GP-hierarchy: 



k 

i9i7P(xfc;x',) = 5^(-A.^. + A,,)7f (x,;xl) 

(1.24) 

j dxk+iiS{xj - Xk+i) - 5{x'j - Xk+i))'yi''~^^\xk,Xk+i;x'i.,Xk+i 



with the correct coupling constant a = Svrao. Notice that if u{t,x) is a solution of the GP equation 
(HHJ), then 

k 

7f ^(xfc;xfc) = Yluit,Xj)u{t,x'j) (1.25) 
i=i 

is a solution of the hierarchy (|1.24l) . To conclude the factorization property ()1.21|) . it remains to 
answer the following three open questions. The first one is to rigorously justify the cutoff modification 
of the Hamiltonian by controlling the clustering of particles. The second one is the uniqueness of 
the solution of the hierarchy in a certain space. The third one is to prove a-priori bounds on the 
density matrices of the Schrodinger equation so that their limits fall into the space needed in the 
uniqueness theorem. Recently the uniqueness problem was solved in Sj. Furthermore, the required a- 
priori estimates were established for certain mean field Hamiltonians without the cutoff modification 
(see Section l3.ll for more details). However, the a-priori estimates and the removal of the cutoff 
modification for the GP scaling remain interesting open problems. 



2 Definition of the Model and the Main Result 

We now define the modified Hamiltonian and state the main result. Throughout the paper we use 
the notation that a ^ 6, if a/b = 0(A^^") for some q > 0. The notation W^'"'^ will stand for the 
standard Sobolev spaces. 



2.1 The Two-Body Problem 

We consider the problem of a single particle in the field of the scaled potential 

Va{x)/2 = (l/2)iVV(iVx). 

Here we assume V{x) to be a smooth, spherically symmetric, and compactly supported potential. 
For < K ^ 1 let and (1 — w'^) be the lowest Neumann eigenvalue and eigenfunction on the ball 
{y : \y\ < k}, 

(-A + iK)(l-u;'^) = e«(l-u;'^) (2.1) 

with normalization ty^dxl = k) = 0. We can extend to be identically zero for \x\ > k so that w'^ 
satisfies 

(-A + iK)(l - w^") = e,(l - w;")xk(x), (2.2) 

where XKix) ■= x{\x\ ^ ^) is the characteristic function of the ball of radius k. We will prove in 
Lemma lA.ll that 

e« = 3aK-3(l + o(l)), 

as long as a ^ K ^ 1. 
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Let ii be a scale to be fixed later on with a <^ ii <^ 1. Let ^{dn) := g{n)dK be a probability 
measure supported on [£i/2,3ii/2] with smooth density g. Define w by 



l-w:= y (1 - w^)n{dK) . (2.3) 

One can check that w satisfies the equation 

{-A + ^Va){l-w)=q{l-w) (2.4) 

where 

/e^l(|x| < K){l-W^{x))fl{dK) 

•= m ^^ M \ • ^^-^^ 

J (1 — w'^{x))iJ,[aK) 

Some properties of the functions w{x) and q{x), which will be used in the rest of the paper, are 
collected in Appendix (see, in particular. Lemma lA.2|) . 

2.2 Removal of triples 

We introduce a new length-scale i with £i <C ^ <C 1. The following procedure excludes configurations 
with more than three particles within a distance i from each other. 
Let h be the exponential cutoff at scale i defined by 



h{x) := e-^^l' , X e A . 

For a configuration x = (xi, . . . x^r) let A/ij(x) be the number of particles other than i and j that 
are within distance £ to either i or j: 

Mijix) := ^ [h{xk - Xi) + h{xk - Xj)] , i j ■ 
Let < e < be fixed. Define 

F(n) :=e-"/^', u>0 

and 

F,,(x) :=F(AA,,(x)) . 

Thus Fij{x.) is exponentially small if jx^ — Xi\ < i or jx^ — Xj\ < ^ for some k ^ The functions 
Fkj, modulo an exponentially small error, forbid particle triples within a distance i and provide a 
strong non-overlap property. This fact and some other important properties of the functions Fij will 
be presented in Appendix IbI 
Introduce the notations 

Xjk = Xjki^) ■= l(ki - Xk\ < h) if ky^j, 
Xjk = Xjk{^) ■= - Xk\ < i) if k j , 
and Xjk = Xjk = if j = A;. Analogously, we will freely use the shorthand notation 

^ij ~ ^a(^Xi -^j ) ? ^'ij ~ ^^Xi -^j ) 1 ^ij ~ ^(^i *^i) ' Q^j ~~ Q^^i *^j) ? 

and similarly for their derivatives. In particular wa = 0, {'Vw)ii = 0. 

We also fix a smooth function 6 £ Cq°{R.) with 6{s) = 1 for s < 1 and 9{s) = for s > 2. For 
some K > we introduce the notation 



^kj 



The constant K will be chosen sufficiently large but independent of at the end of the proof. The 
X-dependence will be omitted from the notation. 
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2.3 The Modified Hamiltonian 



Consider the function w defined in (|2.3j) for the length scale ii. We define 

N 

G,{x) ■.= l-"^w{xi-Xj)F,j{x) , and W :='[[y^ . (2.7) 

The function W is our approximation to the wave function of the ground state. Due to the cutoffs it 
differs slightly from the form mentioned in Section^ It is also slightly different from 

Dyson's construction of the ground state, which is non-symmetric. By Lemma lA. 21 w < do < 1 is 
separated away from 1. We will prove, in Section |S1 (see also Lemma [4.1|) . that there exists a constant 
ci > such that ci < Gj < 1. 

Introduce the notation M^j by 

Mkj:=^{Gj'+G,'). (2.8) 

The modified Hamiltonian H is defined by 

H:=H- Y^i^Vkj - qkj)[l - (1 - Wkj)Mkj] . (2.9) 

Since V and q have compact support of size at most of order £i f Lemma IA.2l in Appendix^, the 
factor {l/2)Vkj — qkj vanishes for any k j unless Xj| < 0{£i). Suppose now \xk—Xj\ < 0{£i). If 
there is no third particle at Xm with m / k,j, such that \xk — Xm\ ^ £/\ log^| then Gj Gk ^ 1 — Wkj 
and Fkj ~ 1 with exponential precision. Thus [1 — (1 — Wkj)Mkj] is exponentially small. This shows 
that the difference between H and H is exponential small unless three or more particles are closer 
than i to each other. Since we will choose i = N~'^^^~'^, for some k > 0, the probability for this to 
happen is, for large A^, very small. 

The reason for this modification will be clear in Section ^ without the subtraction of the three- 
body terms in (|2.9|) we are not able to prove the a priori estimate for (part ii) of Proposition 

EH 

Our methods can be easily generalized to remove only n-body collisions instead of removing 
triples collisions, for any n > 3. More precisely, let /> be a smooth function such that /(x) = 1 
for X < 0, /(x) ~ for x » 1. We define 

and 

4")(x) := FW(Ar,,(x)). 

Then, apart from exponentially small errors, Fj^"\x.) = 1, unless there are at least n — 2 other 
particles in the i'- vicinity of Xj and xj (and Fj^^'^ ~ if this is the case). The modified Hamiltonian 

=H- 5^(114, - q,j)[l - (1 - w,,)Mlf] (2.10) 

with M^f defined through (fTH)) . with Fij replaced by F^^\ Note that now differs from H by 

a term which is exponentially small unless there are n or more particles closer that i = N~'^^^~'^ to 
each other. 
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Since we are dealing with bosons, the Hamiltonian ()2.9|) acts on the Hilbert space 7^®"^ := 
L^(A, dx)®'''^ of functions of 3N variables which are symmetric with respect to any permutation of 
the N particles, i.e., ip E '}{^sN j£ ^^^y ip ^ -^^N j^^^ _ ^ £qj, permutation vr e S^, 
where 

{R^Tp){xi, xn) ■■= ip{x^i, • • • , x^n) ■ (2.11) 
The dynamics is given by the Schrodinger equation 

idti>t = Hipt, or idt'yN = [-H',7Af] , (2-12) 

for a density matrix 7Ar. 



2.4 Weak* Topology for the Kernel of Density Matrices 

Let n := L^{A,dx) be the one particle Hilbert space, and let 'H®''^ be the symmetric subspace of the 
n-fold tensor product H®'^ . We denote by C\ := C^{l-L®^) the space of trace class operators on the 
Hilbert space TL®"^. We will work with families of trace class operators F = {7^"^}n>i with 7("') G C\ 
for all n > 1. 

In this section we consider the operator kernels 7*^"')(x„; x^) as elements of L^(A"' x A"), with the 
norm 

||7(")||2 = j dx„dx'j7W(x„;x'J|2. (2.13) 
For F = {7(")}n>i G e„>i i'(A" x A"), and any fixed z/ > 1 we define the norm 

oo 
n=l 

and we define the Banach space H^^'^ := {F G 0„>iL2(A" X A") : ||F||^(.) < oo}. We also define 
H^^^ := {F G 0„>i L2(A" X A") : lim„_oo J^"||7^"^||2 = 0}. We equip H^^^ with the norm 

||F|| M =supz.«||7W||2. 

n>l 

Similarly to the standard proof of the duality = Cq between scalar valued summable sequences 
and vanishing sequences, one readily checks the duality = (^H^^Y . By the Banach- Alaoglu 

Theorem, the unit ball in H^^ is compact with respect to the weak* topology. 

We shall fix a time T > for the rest of this paper. Denote by C([0, T], H^'^) the space of weak* 
continuous functions from t G [0,T] to Since the space is separable, we can fix a dense 

countable subset in the unit ball of denoted by {Jj}j>i. Define the metric on H^^ by 

oo oo „ 

p(F,r) := 5^2-1^ / dx„dxU(")(x„;x',)[7(")(x„;x'J-7(")(x„;x;)]| . (2.15) 

i=l n=l-' 

Note that the topology induced by p(.,.) and the weak* topology are equivalent on the unit ball 
b'^!!\ We equip C([0, T], i^i''^) with the metric 

p(F,r):= sup p{T{t),f{t)) . (2.16) 

0<i<T 
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2.5 The Main Result 

Consider the rescaled potential Va{x) = {qq / a)'^V {{uq / a)x) (see (|1.7p ) whose scattering length a is 
chosen such that Na = uq. Let £,ii be two lengthscales satisfying 

a<c4<^<l, a£i<^£\ a < £i<^£^/^, £l <^ ai^/\ £^ ^ (2.17) 

(for example £\ = A^~2/3+k £ _ jy-2/5-K for k > small enough). 

Let H be the modified Hamiltonian given in 1)2. 9|) . Let ■0jv^j be a normalized solution of the 
A^-body Schrodinger equation 

idtiPN,t = HiPN,t . (2.18) 

and let 7^^^(xfc;x'^) be its fc-particle marginal densities defined in (|1.9|1 . These density matrices 
satisfy the normalization 

Tr7S = l (2.19) 
for all A^, k and time t. The main result of this paper is the following theorem. 

Theorem 2.1. Assume the unsealed pair potential V is positive, spherically symmetric, smooth and 
compactly supported and assume that q = (87r)^^(||y||i + ||V^||oo) is sufficiently small (of order one). 
Let T > be fixed. Suppose that the normalized initial wave function ipN,o satisfies the energy bounds 

{^Pn,o,HiPn,o) <CiN and {iPn,o, H^^N,o) < C2N\ (2.20) 

Then the sequence has at least one non-trivial limit point and any limit point satisfies the infinite 
Gross-Pitaevskii Hierarchy JTT^j) in a weak sense with a = Svrao- More precisely: 



i) For sufficiently large v, the sequence T^^t = {7^1}fc>i C C([0, T], i?!*^^) is compact w.r.t. the 
topology induced by the p metric. 

ii) Let Too^t = {7,^\}fc>i ^ C{[0,T], B^^) be any limit point ofTjy^t in the p metric. Then Foo,* 
satisfies 

rr(l-A,)(l-A,)7S<C'= (2.21) 
for all i ^ j, i,j = 1, . . . ,k, for all t £ [0, T], and all k> 1. 

Hi) Any limit point Too^t of the sequence Ti\i^t is non-trivial. In particular 

TriS,t = l and Trj!^\ = l. (2.22) 

iv) Let Too,* be a limit point ofTjy^t md assume hr{x) = [3 / Att) r~'^ h{x / r) for any h G Cq°(A) 
with ji\h = 1. Then, for any k > 1 and t G [0, T], the limit 



lim / dxfc+idj;fc+i /ir(4+i - Xk+i)hr'{xk+i - Xj)72Jt(xfc, Xfc+i; x'^, x'^+i) 
r,r'^Oj ^ ^ ' ^ (2.23) 

_. (fc) / . / ^ 

~- "foo,t\^k,Xj,:X.j^,Xj) 

exists in the weak W~^'^{A'' x h^)-sense and defines ^^'''^{xk ^ X j J 5 X j ) QjS Qj distribution of 2k 
variables. 
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v) For any k > 1 and any regular test function j'^^''(xfc; x^) in the Soholev space VF^'°°(A'^ x A'^), 
we have, for any t G [0, T], 

j dxkdxfc/''')(xfe;xfc)7^|j(xfc;x';,) = j dxfcdx';,/''^(xA:; x';,)72Jo(xfc; x'^) 

-'iZl /dxfcdx',j(^)(xfc;x',)(-A,+A;.)7S.(xfc;x',) 

i=i ° (2.24) 

-STrmoV" / / dxfcdxfcdxfc+i /^^(xfc;x'fc) 
,=1^0 J 

X (5(xj -Xfc+i) -5(4 - Xfc+i))7(^+^)(xA:,Xfc+i;Xfc,a;fc+i) . 

Remark 2.2. T/ie main assumption (|2.2Up is satisfied for ipjy^Q = WcpN with sufficiently regular. 
See Lemma \D.1\ in the Appendix for a proof. 

For simplicity, we state the theorem for initial conditions that are pure states. The same result 
holds for initial conditions 7Ar^o G C^{l~i®''^), Tr 7Ar^o = !> with the bounds 

Tr H-iNfi < CiN and Tr H'^-/n,o < C'aA^^ 

(k) 

As pointed out in the introduction, the factorization of the limit point 7^^ ^ depends on the 
uniqueness of the solution to the infinite GP hierarchy 1)2.24(1 . which is an open problem. 

A similar proof holds if we replace the Hamiltonian H l|2.9() with (|2.10() for any fixed n > 3, 
if Q is sufficiently small, depending on n. 

The structure of the rest of the paper is the following. In Section |21 we explain some of the main 
ideas behind the proof of the main theorem. In Section 0] we define an approximation U^^t of Tjv,*- In 
Section IS] we prove energy estimates for H and H^, which rely on some Hardy type and Sobolev type 
inequalities, which are stated and proved in Sectional The energy estimates are turned into a-priori 
estimates for the approximation [/jv,t and for any limit point Uoo,t of UN,t in Section d In Section |H1 
we show how these a-priori estimates can be used to regularize the delta-function appearing in the 
limiting GP-Hierarchy. Finally, in Section IHl we prove Theorem 12.11 Some technical estimates are 
collected in Section^] and in the Appendices. 



3 Idea of the Proof 

We write the solution of the Schrodinger equation in the form ipjy^t = W(l)N,t, where W, defined by 
()2.7|) . is an approximation of the wave function of the ground state of the boson system. We shall 
implement a general idea that consists of two steps. A similar idea in a quite different context lies 
behind the relative entropy method (29) . 

(i) Construct an approximate solution (in our case, W) to the dynamics and factor out this 
approximation from the full solution tp. 

(ii) Derive an effective inequality governing the remaining part (j). In general, one can prove a 
stronger estimate on (p than on ip. This estimate often involves Dirichlet form with respect to 
the approximate solution. 
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In our case, we obtain the following key estimate on (/)jv^f (Corollary I6.2() : 



(k) (k) 

use Uj^ ^ as an approximation of the /c-particle density 7]^ ^ (they will coincide in the weak N ^ 00 



W^\V,VjcPN,t\^ < C (3.1) 

(k) 

for any fixed i 7^ j. We define Uj^ ^ to be roughly the /c-particle density matrix corresponding to 

the wave function (/)jv,t (the precise definition of will be given in Section We are going to 

use as an approximation c 
limit). Then ()3.1() implies that 

Tr(l-Ai)(l-A,)^7);;i<c^ (3.2) 

Let now be a weak limit point of uj^\ which also satisfies ()3.2() . Then, using this bound, 

we can show, by Proposition 18.11 that uj^\{^ki Xj;:x'f^, xj) is well-defined. Moreover we will show 

f Lemma 19. 5p that the limit points of 7^^^ and of coincide, i.e. that 7^\ = uj^\. Therefore 
^^n^ holds true, and the right hand side of ((T^ is well-defined. Note that ((T^ and ((T^ cannot 

(k) 

be proven directly, since we do not have estimates for Tr (1 — Aj)(l — Aj)^]^^; this is why we need 

to introduce the auxiliary densities U^\. 

A more refined calculation is needed to obtain a = Svrao- The idea of this calculation was 
explained in the introduction; it can be made rigorous (Section |^ via the key estimate 1)3. 1(1 and 
certain generalizations of the Hardy and Poincare inequalities (Sectional. To explain (|3.1j) in more 
details, we first review the related mean field models. 



3.1 The mean-field model 

The mean field Hamiltonian is defined by 

N N 

HN = -Y,Aj+N-'Y.V^-\xj-Xk) (3.3) 
j=i j<k 

where we have introduced a scale parameter r and V^'^\x) = t~^V{x/t). If we fix r and take the 
limit N ^ 00 (the mean field limit), then with a product initial wave function Hepp jllj and Spohn 
j27j proved that the one particle density matrix of the solution of the Schrodinger equation ()1.8() 
converges to the Hartree equation 

idtut = -Aut + (yM * \ut\^)ut . (3.4) 

If we take r — > in this equation, we recover the GP equation (|l.lj) but with the so-called mean-field 
interaction constant a = bo = JV instead of o" = Svrao involving oq, the scattering length of V. To 
investigate the simultaneous limit r ^ 0, ^ 00 in the many body problem, we set r = N-f^. If 
< P < 1, then the scattering length of the interaction potential, N~^V^'^\ is much smaller than 
the range of the potential. So we are in the mean-field regime and the one particle density matrix 
of the solution of the Schrodinger equation (|1.8() is expected to converge to the GP equation with 
a = bQ. If /3 = 1, then the Hamiltonian in (|3.3j) recovers the Hamiltonian in the GP scaling (|1.6j) 
and the limit is expected to be the GP equation with a = Svrao. 

In a joint work with A. Elgart we prove that, for product initial data and < /3 < 2/3, the 
density matrices of the Schrodinger equation H1.8|) converge to a solution of the GP hierarchy H1.24() 
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with a = bo- In [3], we also prove that the GP hierarchy has a unique solution in an appropriate 
Sobolev space. Moreover, we show in jS] that for < /3 < 1/2 any limit of the solution to the 
Schrodinger equation with product initial data satisfies the a-priori Sobolev bound. Hence, in this 
case, we prove the uniqueness as well, and therefore we complete the derivation of the GP equation 
with a = bo in a, mean-field scaling. The same result is expected to hold also for 1/2 < /3 < 1, but 
the a-priori bound in this regime is open. 

The outline of the proof in ^ is similar to this paper's. The first step is to prove the estimate 

^\y^yjMt\^ <CN\ (3.5) 

The proof of the convergence to the hierarchy makes use of arguments similar to the ones used in 
Section IS] and Section |21 Since the proof of the estimate 1)3. 5() is instructive, we reproduce it here. 

Proposition 3.1. Suppose tpix) is a smooth normalized function. Then the following inequality 
holds 



(V,//^V) > (l - E / dx|V,V,V(x 

x;/dxivMx^- ^^^^ 



(3.6) 



dxl^fxlP. 



In particular, for t = N ^ with < /3 < 2/3, we have 

N 

/ dx|V,-V^^(x)|2 < C{^, Hl^) + CN^ . (3.7) 

Proof. By expanding the square of the energy operator and dropping some positive diagonal terms, 
we have 

WRnM? dxA,V;(x)A^^(x) 

+ iV E E / dx(-A,V5(x))y«(x, - xM^) (3.8) 

j=l l^m 

1 ^ r 

+ nT.T.J dx(-A,^(x))yW(x, - x^)^(x) . 



j=l i^m ' 

From integration by parts, the second term on the right hand side of the previous equation equals 

N ^ , N 

N 



^EE /dx(-A,V^(x))yW(x,-x™)V^(x) = l^ J] |dx|V,V(x))|VW(x,-x™) 

j=l ij^m j=l iy^m 

2 ^ r 

+ TV Ey dx(V,V5)(x)V^(x)V,yW(x, - X,) . 
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From the Schwarz inequality, the last term is bounded by 



dxVjV'(x)V'(x)Vjy(^)( 



1 



r3/2 
1 



T 



3/2 



dx|VjV'(x)| 

dx|^(x)|2i 



2 1 



vv 



vv 



Xj—X£ 



By Lemma l5.2| the last two terms are bounded below by 



'r3/2 



TV 

■ ^ / dx(|V,VjV(x)|2 + |V,V(x 



)P + lV'(x)P 



The last term in (j!-{.8j) can be bounded in a similar way. Combining these equations, we have proved 
Proposition 13.11 □ 

To understand the restriction r ^ A^^^/s^ consider a single particle in the potential N^'^V^'^\ 
Define the operator P) in by 

t) :=-A^ + N-W^^\x). 

Clearly, (/, f)^/) is given by 



(3.9) 



Notice that for a general smooth function / living on scale one, the last term (/, A^~2(y(T))2j^ 
of order N^'^t^^. This term diverges in the large N limit if r ^ A^^2/3_ ^pj^g middle terms, on the 
other hand, are bounded by CN~'^. Hence for smooth functions living on scale one, (/, t)^/) typically 
diverges. Therefore, functions with (/, P)2/) finite (in particular, the eigenstates of f)) should have 
short scale structure. For such functions, the divergent parts of A^~2(y(T"))2 (^j^j that of A A) are 
cancelled by the middle terms AA^"^^^'^) + N^^V^"^^ A. This cancellation can be understood using 
the following idea of "resolution of singularities" with the help of an approximate solution W. We 
explain the method in the setting of our modified Hamiltonian H. 

3.2 Resolution of singularities 

In this section we explain the idea behind the estimate (see p.23|) ) 



E 



(3.10) 



Since the right side is a constant of motion, this implies the a-priori bound ()3.1() . Recall that H 
contains the singular potential N'^V{Nx). The difficulty of the singular potential was resolved in 
the work of Lieb-Yngvason for the ground state energy problem. The basic idea was to replace 
the singular potential by some flattened out potential and to use the variational principle and the 
positivity of the potential. Our aim, instead, is to estimate the Sobolev norm of (/) in terms of H^. 
Therefore, there is no variational principle and the usage of positivity of the potential is limited. 
Our strategy is to factor out the approximate ground state W. The resulting effective Hamiltonian 
contains a flattened interaction potential, similar to the one in |21]. In this way, we also maintain the 
almost perfect cancellation between the potential and kinetic energy operators, a critical property 
explained in the previous section. 
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The action of the Hamiltonian H on a product function ip = W(j) can be written in a more 
convenient form. Let B be the function 

B:=W-^HW. (3.11) 

Define the operator 

L-=Y,[-^k- 2Vfc(log W)Vk) ■ (3.12) 

k 

The operator L is self-adjoint with respect to VF^, since 

N 



I 



W^^l{L<l,2) =Y1 j W^'(Vfc<^l)(Vfc02) = / VF2(L^i)02 
fc=l 



Then we have 



W-^H{W(l)) = L(j) + B(t), j ^H^ = Y^ j W'^\Vk<t)\^ + j W^B\(j)\' 

k 

and the fohowing identity 



\HW<I)\^ = j W^\L(f)\'^ + j W'^B<j)L(j) + J W'^B(l)L(j) + J \(j)\'^{HWf . 
Last equation can be rewritten in the fohowing more convenient form: 
J \HWcl)\^= J W^lViVjcPl^ J W\V^VjlogW)V,^Vjcj) 

+ 2^y W^B\Vm<P\^ + J W\VraB) (0V™0 + 0V™0) + J B^\cl)\^W' 



(3.13) 



For more details see (|6.15|) - (|6.17|) . If W is the true ground state of H, then i? is a positive constant 
and the middle term in the second line of the last equation vanishes. From the Schwarz inequality, 
we have 

[ I W0p > ^ / VF^|ViVj(^|2 - 2^ / W^\V^VjlogW\\V^^\'^. (3.14) 

The singularities of the ground state wave function are of the form l — Ca/\xi — Xj\ for some constant 
C and for r < \xi — Xj\ <^ 1, where r is the radius of the support of Va {r is of order a). The second 
derivative of the ground state wave function is bounded by Ca/\xi — Xj\^ < C Q/\xi — Xj\'^ in this 
regime, since a < C/K < C'o~^||l^||oo'"^ < Crg. Recall that q = (87r)~^(||l/||i + ||F||oo)- The same 
estimate also holds for \xi — Xj\ < r (see Appendix A). Thus |VjVjlogW^| < C Q/\xi — Xj\'^ . Applying 
the Hardy inequality, the last term can be bounded by 

W^\ViVjlogW\ |V,0p <CgY^j W^lViVjCpl"^ . 

Assuming that g is small enough (of order one) , together with 1)3. 14(1 we get the key estimate 1)3. 10(1 . 

In practice, we do not know the exact ground state function and its approximation will be used. 
When Xi is near xj, the approximate ground state behaves like the ground state of the Neumann 
boundary problem given in (|2.4|) . Thus the singularity of B for |xj — Xj| <C 1 behaves like q{xi — xj) 
where q, defined in (|2.5|) . lives on scale ii » 1/A^. This procedure, roughly speaking, replaces the 
singular potential Va by an effective potential living on a bigger scale li. The price to pay is the two 
middle terms in (|3.13() . Since the singularities of W and q are milder than that of Va, this procedure 
in a sense resolves the singularity of Va- This is the key idea behind the proof of the estimate ((3. If) . 
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4 Construction of the approximate solution 

Recall the definition of the functions and W, 

N 

G^{x) :=l-Y,Hxi-xj)F,,{x) , and := J] . (4.1) 

jfti i=l 

In the next lemma we prove that Gi is separated away from zero. 

Lemma 4.1. There is a positive constant c\ depending only on V such that for sufficiently large N 

ci<Gi<l. (4.2) 

Proof. This can easily be seen because, for fixed i, the sum '^ij-^ij equals zero, if \xi — Xj\ > i 
for all j 7^ i, while, if \xi — Xm\ < ^ for some m ^ i, it is bounded by 

X{\xi - Xm\ <^)Y1 '^ii^v ^ WimFim + Ne'"^'' < ci + Ne'"^" 

where we applied Lemma IB. II So if N is large enough we obtain ()4.2|) with some ci S (cq, 1). □ 

We will also need versions of Gi and W which are independent of some of the variables xi, . . . a; at. 
We define, for i, m 7^ /ci , . . . , A;^ fixed, a version of Fim independent of Xk^ , • • • , x^^ : 

Ft-'-^ := FiMt-'--^ ^ft-'"-^ ■■= E + hn.u] . 

uj^i,m,ki,...,ka 

Moreover, for i ^ ki, . . . ,ka, we define 

(-i{kl-ka) I 1 _ \^ rpiki...ka) 

\ m^i,k\,...,koL 

and 



^(k,...k^) .^^ -Q ^Gf'-''"\ (4.3) 

a 

Thus is independent of xj-i, ■ ■ ■ ,Xk^- When a = 1, we shall drop the label a = 1 and 

denote the function by W^^\ Notice that we have 

Gf ) =G, + W.kFik - ^irniF^S - Fim) ■ 

m,j^k,i 

We shall also use the notations VF^*^! := F-'^ := F^^^'^'^ for any i,m > k, and Gf^^ := 

G\ ■" , for i > k. In Appendix O we shall prove the existence of a constant C, such that, for any 
choice of the indices ki, ..k^, 

C^" < W^'^^-^'^^ /W < G" . 
Since W ^ (by Lemma l4.1j) . we can write the solution of the Schrodinger equation (|2.1(S|) as 

tpN,t = W(t>N,t- (4.4) 
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We define a modified A:-particle marginal distributions by 

Ul^lixi,. . . ,Xk;x[,. . . ,x'f,) := / dxk+i...dxNiW^''\xk+i,...,XN)f 

J _ (4.5) 

X(t>N,t{^l, ,Xk,Xk+l, ■ ■ ■ ,XN)4>N,t{^l^ ■ ■ ■ ■ ■ ■ ,Xn) ■ 

These density matrices will satisfy the requirements of the Main Theorem 12.11 

Clearly ul^\ > as an operator on L^(A, dx)®"^. Furthermore, by (|C.2|) . we have 

■Xk,Xk+l..XN)\ 

< [dx,... dxM l^(x)2|</.(x)p <C'' [ dx|^jv,t(x)P = . 



It is instructive to compare the approximate densities U^l with the true marginal densities 7^^^, 
which are given by 

.... .4) := /d...« . . . d« wi., .... .... 

X 4>N,t{^l, ,Xk,Xk+l, ■ ■ ■ ,XN)4>N,t{^'i' ■ ■ ■ • • • ^^n) ■ 

A simple computation shows that jj^^ and satisfy the relation 1)1.19^ . 

5 Hardy type and Sobolev type inequalities 

We need the following finite volume version of the usual Hardy inequality 

I %<c[ |V/|2. (5.1) 

J«i \xr Jr3 

Lemma 5.1 (Hardy Type Inequality). Let <U{x) < c|a;|~^ on B := {\x\ < £} C M^, and let 

if) -.= 131-^ [ f 
Jb 

be the average of f on B (\B\ = {A/3)tt£^ is the volume of B). Then 

{U\f\')<C{\Vf\') + C{U){\ff). (5.2) 

Proof Let / := (/), then 

{U\f\')<2{U\f-m + 2f-^{U). 

By the usual Hardy inequality 

- /T) < c(|v/|2) + cr^d/ - . 

By the Neumann spectral gap we have 

r2(|/-/|2)<c(|v/|2). 

The lemma now follows from 



\fT = \{f)\' < {\f\' 
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□ 



We also need some well-known inequalities, collected in the following lemma. 
Lemma 5.2 (Sobolev Type Inequalities). The following two inequalities are standard, 
i) Suppose U G L3/2(A)^ U>0, and e W'^''^{A). Then 



m''U<C\\U\\Ls;, 



(A) 



(5.3) 



ii) Suppose U S L^{A). Then, considering U{x — y) as an operator on the Hilhert space L^(A, drE)c 
L^(A, dy) we have the operator inequality 



U{x-y)<C\\U\\L. (1-A,.)(1- Ay) . 
Proof, i) Applying the Holder inequality we find 



(5.4) 



1/3 



By the Sobolev inequality we have 



1/3 
|6 \ < 



which implies part i). The proof of part ii) can be found in [Hj. 
Let 



\{x) 



X{\x\ < 3£i/2} 



and a{x) :- 



X{\x\ < 3h/2} 



X e A , 



□ 



(5.5) 



and, as usual, Xij := A(xj — xj), CTij := (j{xi — xj). We note the following inequalities, whose proof 
is given in Appendix IXI 



w < Xj "f^ ^ CaiiX, \Vw\ < Ca 



-1 



\Vw\ < CaX, 



\V^w\ < CqX, \V^w\ < Caa, \Vw\^ < Caa, < Ca^X, ^^'^^ 

where i)= (87r)-i(||y||oo + llV^Hi). 

In the sequel we will use the convention that integrals without specified measure and domain will 
always refer to Lebesgue integration on A^: 



. .. dx . 



We will use the Hardy-type and the Sobolev type inequalities in the following form: 

Lemma 5.3. iei (/>(x) G L^(A^). For any q > and for fixed k, j indices and fixed xi, ..,Xk, ..,xn 
we have 



I dxkW^F^^XkM^ <C, I dxkW^F'^^%,\Vk^\^ + C,£ir^ j dxkW^F'^^%j\(t>\^ 



for sufficiently large N. In particular 



kj k 



\4>\' 



(5.7) 



(5. 
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by \B.^) . Moreover we have 



J W^akM^<C{logN) J [^|V,V,</.|2 + iV^|V,-</.|2 + iV 



2| A|2 



(5.9) 



kj \ jk 

The same estimates hold if we set W = 1 in all integrals. 
Proof. We define 



F 



(k) (k) 

Note that is independent of the variable x^, and that, trivially, Fkj ^ Fj . Moreover we have 



Xkj Fl- = Xkj e 
> Xkj e 



-1^ ' Tlm^k,]W^j-^"i) + h{Xk-Xm)) 



(5.10) 



> Xk, [FfY. 

because 1 + e < 4. Here we used that ((x^ — x^)^ + f'Y^'^ < {{xj — Xmf' + ^^)^/^ + \xj — Xk\. 
Using ()C.1|) to replace W"^ by [W^'^'^^]^, which is independent of x^, applying 1)5.2(1 for Xjt, and finally 
changing the measure back, we find: 

<C [ dxi...d^k---dxN [W^(^)]2(Ff V 



dxfc 



Xkj 



<c I w^F^'^Yxkj \VkcP\'' + hr'\cp\ 



-3U|2 



The estimate ()5.7() follows now from (|5.1U() . Similarly the estimate ()5.9() follows from ()5.4() . after 
replacing W by W^'''^^ It is clear that setting W = 1 would not alter the validity of the proof. □ 



6 Energy Estimate 

We begin by computing the action of the Hamiltonian H defined in ()1.6() on the wave function W, 
defined in Section [T^ For any fixed index k we have 



Here, and in the rest of the paper, we use the summation convention: all unspecified indices are 
summed up; in this case the j and i indices on the right hand side. 
Direct computation shows that for any fixed k 

-^kGk = {Aw)kjFkj + ^k , 

where 

^k ■= '^C^w)kjVkFkj + WkjAkFkj (6.1) 

and, if 7^ j are fixed, then 

-AkGj = {Aw)kjFkj + ^kj 
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with 

:= 2{Vw)kjVkFkj + WijAkFij . 

From it follows that 

Aw = -(1/2)K(1 -w) + {l- w)q . 

Then we have 

W-^HW = {{l/2)Vkj - qkj)[l - (1 - Wkj)Mk,\ + quj + ^ (6.2) 

with 



and, for any (see 12? 



:=^(G-'+Gl-'). (6.4) 
The action of the modified Hamiltonian H2.9() on the wave function W is then simply given by 

W-^HW = qkj + ^. 

It follows that the action of H on the product W(j) is given by 

W-^HW^ = L^ + B^ (6.5) 

where we defined the operator 

L := ( - Afc - 2Vfc(log W)Vk) , (6.6) 

k 

and the function 

B := VF-i W = qkj + . (6.7) 
Note that the operator L has the self-adjointness property that (with the summation convention) 

W^iPiLcj)) = j W^{Vki>){^k<i}) = j W^{Li^)4> . (6.8) 

Next we compute the last two terms in (see (|6.3j) ). For fixed i,k we have 
VfcGj = {Vw)ikFik - iVw)k£Fki5ki - WiiVkFi£ . 

Hence we obtain 

(VfcGi)^ = {Vw)fkFfk + {Vw)ke{Vw)ksFkeFksSki - 2{Vw)ikWiiFikV kFm 
+ '^{^w)keWksFkeV kFksSki + wuWigV kF^V kFis . 

Using Lemma iB.ll we note that the second term on the r.h.s. is exponentially small unless i = s. 
Analogously the third term is exponentially small unless k = i, while the fourth and fifth terms 
are exponentially small unless i = s. Thus summing over i,k on the left and the right side of the 
equation, we obtain (changing the name of the indices in the second and fourth term) 



jVkG.f _ iVw)lFl {Vw)ikW^kFik{VkF^k - V^F^k) . wU^kF^if ( ,-ct-^\ (aa\ 
4G2 2{l- WiuF.,k? 2G] + 4G2 + I' J' ^^'^^ 
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Here we also used that, on the support of {'Vw)ik, we have Gi = l — WikFik + 0{e ^1^''). Analogously, 
for i ^ j we have 



^GiGj 2{l -w,kFik? 2G,G, AGiG, 

(6.10) 

where we are summing over all k and all i ^ j on the left and the right side of the equation. We see 
that there are some cancellations between 1)6. 9|) and (|6.10() . It follows that 

n = n + 0(e"^~'^ , (6.11) 

where 

f^ = ^ + ^ + r (6.12) 

and 

' (yw)ikWjeFikiVkFje - ViFji) wuWjs^ kFuV kFjs 



r — 



2GiGj AGiGj 



Here Sij = —1 \i i = j and Sij = 1 \i i ^ j (and sum over all indices is understood). A more careful 
analysis also shows that, analogously to 1)6. 11() . for any fixed m, 

Vm^ = V„Sl + 0[e-^^~"^ . (6.13) 

This fact will be used in the proof of the next proposition, which is the key energy estimate. 

Proposition 6.1. Assume a = a^N^^ , a <^ £i <C ^ <C 1, aii <^ i^, a <^ i'^^'^ and o <C (for example 
£^ = Ar-2/3+« and £ = N-^'^-^, for k > sufficiently small). Put g = (87r)-i(||y||oo + \\V\\i). 

i) Then we have 

{W^,HW^) = j W^{HW(t))>{l-o{l)) j Ty2^|Vfc(/)|2-o(iV) / W^\(t)\^ , 

k 

for N ^ oo. 

ii ) Moreover, if g is sufficiently small, there is C > such that 

J |W</.|2>(C-o(l)) J W^2^|V,V,(A|2-o(l) (n I iy2^|V,<A|2 + Ar2y W^\cpA 

ij \ i / 

as N ^ oo. 

This Proposition and the conservation of H and along the solution of the Schrodinger equation 
(|2.18|) immediately implies the following 

Corollary 6.2. Let '(/'7V,t(x) = W{x)(j)]\f^t{^) be a symmetric (with respect to permutations of the 
particles) wave function solving the Schrodinger equation i2.18\) . with initial data ipN,o satisfying the 
hounds 

{^Pn,o, Hi^Nfl) < CiN , and {iPn,o, ^ViV,o) < C2N^ . 

Suppose moreover that the assumptions of Provosition are satisfied. Then for any time t and for 
any fixed indices i ^ j, we have 

J W^\Vi(j)N,t\^ <C, J VF2|ViVj0iv,t|' < C (6.14) 
for some constant C, independent oft and of the indices i,j. 
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Proof of Proposition \6.1l From (|6.5() , (|6.8() we have (recall the summation convention) 

j w^{Hw^) = j w^\Vk<l)\^ + j w'^B\(i)\^> j w^\Vk^\'^-j w^inM"^ . 

Here we used that B = q^j + 0, > 0, because qkj > 0. Part i) thus follows immediately from H6.11() 
and Lemma 16.31 below. 

As for part ii) we note that, from (|6.5|) . 

j \HW(j)\'^ = j W^\L(j)\'^ + j W^B^L(j) + j W^B(j)L^ + j \(l)\^{HWf . (6.15) 

Using 1)6. 8() . we can rewrite the second and third term in the last equation as 

f W'^B^L(^+ f W'^B(f)L^ = 2 f W^B\V^^\'^+ f W^VmB) (0V„0 + ^V^,^) 

(6.16) 

>2 / W^n\V^6\^-2 / t^2|v R|U||v^,(*|. 



Here we used again the positivity of q. From Eqs. ()6.11|) and 1)6. 13() . we find 



W^B(I)L4>+ / W^B(l)Lct> 



> - '^ j W^M\Vm(t>?-2 j W^\V^qkM\^ra(f>\ 

-2 j W^\vJlM\Vm(t^\+0(e-''~') {^j VF2|v^0|2 + iv j wM 
Analogously, the first term on the r.h.s. of ()6.15)) can be written as 
j W^\L(t)\'^ = j W^Vj{L^)VjCp 

= j W^\ViVj(j)\^ -2 j W^{ViVjlogW)Vi^Vj(l) . 



Next we note that, for any fixed i.j, 



V. V, logW = lYl log Gt = lYl 



2^ ' ^ ^ 2^ GJ 



(6.17) 



The second term gives a positive contribution. Hence we have 

Thus, using > ci > 0, we obtain 

j \HW(t>\''> j VF2|v,V,(/.|2 - (1/ci) j W^\ViVjGf\Y^MV,(t^\ 

-2 j W''mV„M^-2 j W^\VmqkjU\V^m<t^\ (6.18) 
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The first term on the r.h.s. is the positive contribution we are interested in, all the other terms are 
errors terms which we estimate separately in Lemmas 16.31 (note the remark after Lemma [6.3p . 16.41 16.51 
and 16.61 The proposition follows then from last equation and from the results of these lemmas. □ 

In the rest of this section we use the following notations for brevity 

k 

kj k 

where the dependence on N is omitted. 

Lemma 6.3. Assume a ^ li <^ I <^ 1 and aii <C i^. Then 



I 



w'^\n\\^\'^ < o(i)f(0) 



as N ^ oo. 

Remark. Replacing (j) by 'Vm<P, and summing over m, it also follows from Lemma 16.31 that 

Proof. Using Gj > ci > 0, we have from (|6.12j) (using the summation convention) 
W^\nM^<C J W^[\{Vw)kj\\VkFkj\+Wij\AkFij\ 

+ C j W"^ \{Vw)ik\WjrFik\VkFjr\ +WimWjr\VkFim\\VkFjr\ jc/ip 

<c j w\ae~' + aeie-^)XkjFlj^\<p\^ , 

where we used the estimates Wjr < Xjr, Wij < CaiiXij and |V^«^j| < CaXij from 1)5.6(1 . and we also 
used the formulas ()B.3|) . ()B.6|) to sum up indices. Using £i < i and ()5.8|) . we obtain 



0-3| i,|2 



□ 



Lemma 6.4. Assume a <^ ii <^ £ and aii <C Then there is a constant C such that 
j ^' 5^ |ViVjG™||Vi</>||Vj,/.| < Cgj VF2^|V,V,-0|2 + o(1)^(<^) 

i,j,m 

for N ^oo. Recall that g = {87r)-'^{\\V\\^ + ||F||i). 
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Proof. A direct computation of 'Vi'VjGm and an application of Schwarz inequality (to separate | Vj(; 
and |Vj(/)|) show that (with the summation convention) 



/ 



W^\ViVjGJ\V^(t>\\VjC 



<Aj W\V^w),,\F,,\V,<l)\'' + C j W\Vw)uj\\V^Fkj\\V,(t>\'' (6.19) 
+ C/ l^^|(V.).,||V.F.,||V,^|^ + / T^VJV.V,F..||V. 



Applying \Vw\ < CaX (see H5.6() l. (|B.3|) . and (jS.Hp . the second term is bounded by 

j W\Vw)k,\\V^Fkj\\Vi<t)\'' < Car' J W^XkjFl^.^lV.cPf < o(l)^(0) (6.20) 

if a£i <C and a <^ i. The third term on the r.h.s. of (|6.19|) can be bounded analogously. As for 
the fourth term in (|6.19|1 . we note that, using (|B.3() and w < CaiiX (by (|5.6j) ). it is bounded by 

I W^Wkm\ViVjFkn,\\Vict)\^ <Cae,r^ I W^\k^Fl^\Vicl)\^ 

and the same estimate holds as in ()6.2U() since aiii"'^ < a£~'. 

Finally we consider the first term on the r.h.s. of (|6.19|) . Using iV^wj < CgX from 1)5. 6|1 and the 
estimate (|5.7j) we obtain 

J W\V^w)ij\F^J\V^(|)\^ <Cg J W^X^JFiJ\V^cj)\^ 

<CQj W%,Fl/^\ViVj<p\^ + CQhr' j W^XijFl/^\Vi<p\^ 
where we also used ()B.2|) . This completes the proof because li^'^N^' = o(l), as ^ oo. □ 

3/2 

Lemma 6.5. Assume a <^ ii <^ i and a <^ . Then 

m,k,j 

for N ^ oo. 

Proof. By part iii) of Lemma IA.2I we have (with the summation convention) 

J J I I 

Next we apply a weighted Schwarz inequality: 



<Caaq'' I ^2 X(|xr^ Xj\< lh) ^^^^ 
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In the first term we use (|5.4() from Lemma 15.21 in the second one we estimate xil^m — Xj\ < < 
Cii\xm — and apply the usual Hardy inequality (|5.1() . In both cases we first have to change 

W to W^^^ by (frmi then back again to W to make the weight function W independent of the xj 
variable. The result is 

Calf I W\{\xm - x,\ < lhM\Vm<t>\ < C + ^ H^k) . 

We choose a = ^^"^ . Then 

Ca^r'/ W\{\xm-Xj\ < \h)\^\\Vm<\A<o(X)T{<ii) 

provided a -^l^ . □ 
Lemma 6.6. Assume a < < £ < 1, a<^l^ and ali < t'^. Then 

as N ^ oo. 

In order to prove this lemma we need the following estimates (we recall the definition of O^rn from 
Lemma 6.7. Let a <C -^i <C ^ ^ 1 and a <^ l"^ , then for any g > 



j W^ak^FlJ(t>\\Vm<t)\ < o{Ni)T{(t>) , (6.21) 

I W^\kjFl.ekmm\Vrr,(t>\ < o(l).F(0) . (6.22) 



Proof. With a Schwarz inequality we obtain (using the summation convention) 
W^akmF',J^\\Vm(p\ < [ W^akmFUa\cl)\^ + a-^\Vr 



<aj W^OkM'' + a-^a-^ J W^XkmF^jV 



where we used that akm < a ^^km- Now we apply (|5.9|) and (|5.7|) from Lemma [5^ to the r.h.s. of 
the last equation. We find 

J W^akmF^J(^\\Vm(l>\ < Cg {alogN + a-^a-^ + a-^ii-^) J^{cl>) . 

Eq. 1)6.211) follows choosing a = for sufficiently small 5 > 0. 

As for ()6.22p we have, from ()5.7p . 

j W^\kjF'^.ekra\m^<t^\ < j W^\k,Fl.ekm{a-^\yra(f>? + a\<l)\'') 

< I W^XkjF'J.^Okm 

X (a- VjV^^P + a-Hir^\Vm(l>\'^ + a|Vj-(/.|2 + aiir^\(P\'^) . 
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Applying ()5.3() for the second and third terms, and 1)5. 4() for the last term, we get 

< C (a-i + a-^lir^i'^\ log£|2 + af\ \og^\^ + ahr^l^\ log^p) ^((/)) . 
Choosing for example a = we find H6.22() . 

Proof of Lemma \6.(A By the definition of we find (using the summation convention) 



□ 



w^\Vmmmm(t>\ < I 



+ / w 



'2Gk 



|0||V™(/>| 



jk 



|V^(/.| + / W 



(6.23) 



IV. 



We begin considering the first term on the r.h.s. of ()6.23|1 . 

/ <C J (|V„l^fc| + \nk\\VmGk\) mVm^l^l . (6.24) 

From the definition of 0,^ (see 1)6. ip ). using Lemma lB.2l to bound the derivatives of F^j, we find 

J W^\Vrr.nkM\VrrM<Ci-' J W\V^wUj\F^.\^\\Vm(l)\ 

+ cr^ j w\vw)kj\Fl'^{e^k + emjm\^'m(t^\ (6-25) 
+ cr' j w^wkjFlj^emk + OmMl^'m^l + o(£^-i-^).F(0) . 

We use that \{V'^w)mj\ < Caamj, \{'^'w)kj\ ^ CaX^j and Wkj < CaiiX^j (see ()5.6() ). and applying 
Lemma 16.71 we obtain 

J W^\VAM\Vm^\ = o{l)T{cP) . (6.26) 
For the other term in 1)6. 24^ . using the definition of 1)2. 7|1 . we have 

J W^\Qk\ \VmGk\ 101 \Vm(l)\ 

'w)km\ Fj^jyi \^ kFkj\ 

+ \(yw)mj\ WmjFmn \^mFmj\ + |(Vtt;)fcm| W^jFj^m \^kFkj\ 
Wkn \^mFkn\ \^kFkj\ WkjUJ^n 

Taking as example the first term on the r.h.s. of the last equation, we note that, unless j = n, 
the term is actually exponentially small because, by Lemma IB. 11 ||xmjXmn-^mn||oo 

< exp(-cr^) if 

j 7^ n. Applying similar arguments to all other terms as well, and applying Lemma lB.2l to bound 



26 



the derivatives of F, we obtain 



W'\nk\\yraGk 



<C j iy2|^||V„(/.|(rl \{Vw)mn?Fmn + [i'^ \{V w)ki\Fkj + T^W^Ekj) 



o(i)r(4,) 



-cl- 



+ 0(f 



following the estimate of the term ()6.25|) and using the bounds (Vtt;)^ < Caa and w'^ < w from H5.6|l . 
So, together with (|6.26|) we find the following estimate for the first term of (|6.23jl 



'2Gk 



(6.27) 



We consider next the second term on the r.h.s. of ()6.23|1 . Clearly we have 



\VmCt>\ <C J W^ilV^Qkjl + \nkj\\VmGj\) \cj)\\Vr 

Comparing Q^j with 0,^, it is clear that it only remains to control the two terms 

J W^W^j\AkF,j\\VmGjM\^m'P\ 

and 



/ 



W^\Vm{WijAkF,j)M\Vr 



(6.28) 
(6.29) 



We begin with (|6.28p . The summation over k is performed by ()B.3|) . We have 

W^W^j\AkFij\\VmGjM\Vrr,cl)\ < C j (wijFll\V w) j.m\Fjr 



1 /2 ± z 

+ WimF^J^ \{Vw)mn\Fmn + WijFl!''Wjn\VmFjn\] |0||Vm</'| ■ 

By Lemma IB. II we know that, for example the first term is exponentially small unless i = m. 
Similarly, apart from exponentially small contributions, in the other two terms we can consider only 
the case i = n. Hence we find 



V2, 



W^Wij\AkFij\\VrnGjM\^m^\ < J W^WmjF^^^\{Vw)mj\ F^j |</.| \V ^ct>\ 

+ Cl-^ j W'w%F}l'\VmF,n\ \^\ \Vm^\ (^.30) 

+ O (exp(-cr")) £{(!)) . 

On the right hand side we estimate \Vw\ < GaX hy 1)5. 6() . In the first term, using F < 1, we have 

Cr^ J W^WmjF^^\Vw)mj\ Fmj \cl)\ |V^0| < Cai'^ J W^Xmj Fmj \cl)\ |V„(/>| (6.31) 
that can be easily estimated by o{l)J^{(p), using the assumption a <^ £^ and a <^ ii <^ i. 
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In the second term of H6.30() we estimate < Ca^X. By Lemma IB. 21 we find 



where we also applied Lemma 16.71 
We consider next ()6.29() . We have 



W'\VmiwijAkFij)M\V^<P\ < / t^^|(Vw;)„,-||AfcF„,-||<^||V„</)| 



+ J W^WijlVmAkFijMl'^r, 



(6.32) 



After summation over k using (|B.3() . the first term can be treated as in (|6.31|) . 

As for the second term in 1)6. 32() . we use that w < CaiiX and Lemma lB.21 We find 



which can be estimated by o{l)T{cj)) by Lemma \^7l\ 

Finally we control the third term on the r.h.s. of (|6.23|) . After applying Lemma IB. II we find, for 
any fixed m, 



v^r <c 

+ 
+ 
+ 
+ 



jm\ Fij^ 



jm\ 



\{V'w)ik\ Wjn |VmFi.,, , ^ , ^ 

\iVw)f^\ WjnFfm \^mFjn\ + \{^w)lj\ WjnF^^ |VfcFj, 
\{V'w)ik\ WikWjnFik \^kFjn\ \'^mFik \ + \ {^w)ik\ w'j^F, 
\{^w)ik\ \{yw)jm\ WjmFikFjm \^ kFjm\ + |(Vw)mi| W. 



W)ik\ Wjn |VmVfcFjn| Fik 

"in| + |(Vw;)L ■ 



- KViyjfc^l Wjni'km IVfci^jnl 
IV^Fifcl + \{Vw)ik\ W^jnFik l^kFjnl \V„rFjn\ 
)m \^kFjm\ + |(Vw)mi| Wjn IVfeF^il |VfcFj„| 
zi\\^kFjn\ + WinWjs\y k'^ mFin\ |VfcFjs| 



(6.33) 



lib'- Lts:^ Jilt I ' h:-^ JTYll mi I "^jn I ^ fe-* mi \ \ 

+ \{^w)mi\ WmiWjnFmi | VfcFmj 1 1 VfcFj„ | + WinWjs\V k^ mFin\ |VfcFj, 

+ WinWjs\ykFin\ |V, 



' mFin 



\^kFjs\ 



Using the estimates (|B.3|1 . ()B.6|) and (|B.7|1 from Lemma [6.21 and the bounds ()5.6|) . we obtain 



< C 



ar^GmkFlll + (ar^ + ahr' + a'r'')XikFll\Q.,m + 



-,1/2, 



ik'^ ik \"im T "km) 



+ |(Vtf)ifc| \iVw)jm\ Fik |VfcFj„| + \iVw)ik\ W%Fik \VmFjn\ IVfcFj 



(6.34) 



The terms on the first line of the last equation can be bounded using directly Lemma 16.71 Consider 
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next the first term on the second Une. Its contribution, after two Schwarz inequahties is given by 

W\Vw)ik\\{yw)jrn\Fik\ykFjrnM\ym(t>\ 



< / W^a\{Vw),k\ 



XikFikXjm\'^ kF^ 



]m\ 



+ j W^a-\yw)jrn?y^\(t)\'^ +^~^\V^4>\'^ \ xikFikXjn^\ykFJrn\ 

<Cal3ar' J W^aikFik\cP\^ + Ca(3~H~^ J W^XikFik\VmCp\^ 

+ Ca-^^ar' J W^aj^Fj/^\ct>\^+Ca-'r'e-^ J W^\j,^f]I^\V m<P? 
< C [af5ar^ log N + al3-^r^{l + iiC^) + a"S«^"^ log 

Optimizing the choice of a, j3 and 7, we find that this term is of order o{\)J^{ifi) (using that a <C 
and ah < l^). 

As for the last term on the r.h.s. of (|6.84|) we use that < Xjn- By Lemma lB.21 and Lemma 
15.31 we get 

j W\Vw)ik\ w]^Fik \VmFjn\ \ykFjn\ |V„^(/>| 

< Car^ j W^\ikFik{a\^? +a-^\ym(t^\'') 

< C{a^/^r^ + a3/2^ir).F(</.) = o(l).F((/.). 

From 1)6. 34() . we find 

' 101 |V^</>| = . 



/ 



Vrr,T 



This completes the proof of Lemma 16.61 



□ 



7 A-priori Estimates on Un^i 

Recall that TC = L^(A, dx) is the one particle Hilbert space, and that 7^*^" denotes the n-fold tensor 
product "H®". We denote by £^ := and by /C„ := /C(W^") the space of trace class operators 

and the space of compact operators on the Hilbert space 'H'^'^ respectively. We equip these spaces 
with the trace norm, || • ||i := Tr| • |, and with the operator norm, || • ||, respectively. It is a well-known 
fact that /C* = that is £^ is the dual Banach space of /C„. 

We define Sobolev-type norms on trace class operators. For 7^") G C^{7i^^), we define the norm 



l7^"^llw„ 



Tr|5i7«5i| ifn = l 

Tr|5i527(")5i52| ifn>2 



where Sj = (1 - Aj)^/"^, and the space W„ = {7^") G : h^^^Hw^ < 00}. 
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Note that Wn is the dual of the space 

_ r {SiK^^^Si : G /Ci} if n = 1 

" \ {^i^sKW^sSi : G /C„} if n > 2 

equipped with the norms 

||r«|U, = ||5f ^TW^fi II and ||T(")|U„ := \\S{^S2^T'''''>S^^S^^\\ if n > 2. 

Lemma 7.1. Suppose Ujy^t = {U^\}^^-y is defined as in \4-5^ , and that the assumptions of Theo- 
rem \2.1\ are satisfied. Then there is v > \ large enough such that \\UN,t\\jj('^) ^ 1- Moreover there 
exists a constant C > such that 

\\u^Nl\\y^,<c' (7.1) 

for any t G [0, T] and k < N . 

Let Uoo,t be a limit point ofUN^t in C([0, T],H^^'') with respect to the metric p (defined in \2.1b]) ). 
Then we also have ||?7oo til rrC") ^ 1 o,nd there is a version ofUoot such that 

rr(l - Ai)C/^Jt < C and Tr{l - Ai){l - Aj)U^^l < C'' if k>2 (7.2) 
for every k, every t G [0, T], and every i,j = 1, . . . ,k with i ^ j. 

Remark 7.2. This lemma does not yet prove the compactness of the sequence Uj^^t £ C{^iT]:H'^-'^): 
to this end we still need the equicontinuity ofUN,t, which will he proven in Section \y.2l Here we prove 
that, if a limit Uoo,t G C{[0,T], H^^) exists, then it satisfies j?.'^ . 

Proof. We have 

= / dx,dx',|[/g(x,;x',)|2 < C^TrUi;} < C^K 

(k) k 

Here we used that Uf^[ is the kernel of a positive operator with trace bounded by Cq. this implies 

(k) h 

that also the operator norm of {7j^j is bounded by Cg. Choosing v > 2Co we immediately see that 
II^A^.i llji^C") ^ 1- The bound ||C/oo,t||j|^(i') < 1 follows because the norm can only drop when a weak 
limit is taken. 

From Corollary 16.21 and (|C.2|) we have 

Tr|5iC/^']5i| = j |T^«|2 (|Vi<A7v,t|' + \(t^N,t?) < C J + \<pN,t\^) < C 

and, for any k >2, 

TT\SiS2Ui{lSiS2\ = J \W^''^\ 



|ViV2(A7V,t|^ + |Vi07V,t|^ + \y2(t>N,tf + \(l)N,tf 



(7.3) 



for some constant C independent of k. This proves (|7.1|) . 

For each fixed t, we can assume that UN,t converges to C/oo,t in the weak* topology of H^'' 
(otherwise we choose an appropriate subsequence). This follows because UN,t e B^_^^ (the unit bah 
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of H^^) and because, on B^^\ the topology induced by the metric p is equivalent to the weak* 
topology. This in particular implies that, for every fixed k, converges to U^^^ in the weak 
topology of L^(A^ X A^). Note that U^^^ is symmetric w.r.t. permutations of the k particles, 

ik) 

because ^ is symmetric, and because the symmetry is clearly preserved in the limiting process. 
By passing to a subsequence and using the Alaoglu-Banach Theorem, we can assume that U^^ ^ also 
converges in the weak* topology of Wfc, and let Uoo,t be the limit. Testing these two limits against 
operators with smooth kernels, one easily sees that = as elements of L'^{A^ x A'^), hence 

ik) (k) f/^l 

Uj^^ converges to ^ in both topologies. Then the estimate in 1)7. 2(1 for the version of given 

~ (k) 

by follow from (|7.1() . from the permutation symmetry of U^\, and from the fact that the norm 



does not increase under weak* limit. 

The next lemma will be used in Section |2l to prove part iii) of Theorem 12.11 
Lemma 7.3. Assume that a <^ £i <ti i <ti 1- Then, for any fixed k and t, 

i) We have 



□ 



Tr[U 



\^N,t 'N,t J 







as N 



oo. 



ii) Assume that Uoo,t = {U^\}k>i G C{[0,T], H^^) is a limit point of Un^i = {f/}f|}fcLi with 
respect to the metric p. Then 



TrU 



(1) 



oo,t 



1 



and 



TrU 



(2) 



oo,t 



1. 



(7.4) 



Proof, i) We have (with (p = 0Ar^t(x)), 



Tr 



\ N,t 



(k) 
TNt 



(^2 _ 



j=l j=k+l m=k+l j=k+l j=m+l 



and hence, using HC.6|1 . 



j=l m=k+l 
N 

j<k m=l 



in,j>k r<k ' 

So applying Lemma 15.31 for the first term and part (ii) of Lemma 15.21 for the second term (with 
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u = 

^mr)) we find 




< Ck{a£i + aejr^ + he^i log^p + ej\iogef) + o{e^~') = o(i) 

for ^ cxD using (|6.14|) . 

ii) From part i) we have, for any fixed k and t, Tr — > 1 as A'^ — > oo using the normahzation 
H2.19() . By Lemma [7.11 the sequence in Wi is compact w.r.t. the weak* topology of Wi. By 

passing to a subsequence we can assume converges to in the weak* topology of Wi. In 
particular, 

TV U^^\ = Tr S^^SiU^^\Si = lim Tr S^'^SiU^j}\Si = lim TV U^^\ = 1 

' ' N^oo ' N—*oo ' 

since the operator S'j"^ = (1 — Ai)^^ is compact on the finite periodic box A. The proof of H7.4|) for 

(2) 

U^\- is similar. □ 



8 Approximation of the Delta Function 

We consider the sequence of density matrices UN^t = {U^N\}k>i- By viewing their kernels as distri- 
butions, we will prove in Section IHUl that UN,t is compact in C{[0,T], H^^'^). Assuming this property 
for the moment, we denote a limit point by Uoo,t- Since for each fixed t the kernel Uoo,t is defined 
only as a distribution, the restriction of uj:^~l^^ on the diagonal a^^+i = i.e., 

^^c£j^^(xA:,a;fc+i;x'fc,a;fc+i) (8.1) 

has no a-priori meaning. On the other hand, based upon Lemma l7. II we can view Ui\i,t and {7oo,t (or, 

more precisely, the families of operators with kernels given by C/^) and U^\) as elements of W^'^^ 
for every t £ [0, T]. 

We shall show in the next proposition that since Uoo,t £ W^*^^, the diagonal element 1)8. If) is well- 
defined. For the rest of this section, we shall assume that is the kernel of a density matrix 
with Tr C/C^+i) < oo. 

Proposition 8.1. Suppose (5/3(x) is a radially symmetric function, withO < df3{x) < CI3~^x{\x\ < /?) 
and f 6i3{x)dx = 1 (for example Sf3{x) = f3~^h{x//3), for a radially symmetric probability density 
h{x) supported in {x : |x| < 1}^. Then, for any 

j{k) ^ p^i.oo^^fc X ^fc^) and for any smooth function 
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(8.2) 



[/{'=+i)(xfc_(_i; x^_(_-i^) corresponding to a (k + l)-particle density matrix, we have, for any fixed j < k, 

X {^Pii^'k+l - Xk+l)S/32iXj - Xk+l) - S{x'k+i - Xk+l)5{Xj - Xk+l)) 

< C[\\J\U + ||V,-J||oo] iPi + v%) Tr\SjSk+iU^^+'^S,Sk+i\ . 
Remark 8.2. From H8. using a standard approximation argument, it also follows that 

/ dxfcdx';,dxfc+id4+i J('^)(xfc; x'fc) C/(^+^)(xfc 
X ((^ftl^+i - 2;fc+i)(5/32(xj - xk+i) - (5/3j(xfc+i - xk+i)6f^'^{xj - Xk+i)) 

<C7[ll^lloo + ||V,-J||oo] + + TrlSjSk+iU^'^+^^SjSk+il 

for any U^^^^^ for which the r.h.s. is finite. 

Proof. It is enough to prove for C/C'+i) of the form U^^^'^\y.k+i,^k+i) = /(xfc+i)7(x'fc+i) (in 
this section we use the notation x^+i = (xi, . . . Xk+i) and x'^^-^ = {x'^, . . . , Then, for a general 

density matrix U^^^^\ the proposition fohows by considering the spectral decomposition: 



U^''^^\xk+i;x'^.^i) = ^c„/„(xfc+i)/„(x^ 



k+l) 



with c„ > for all n and ||/n||L2 = 1, and using the fact that Tr [/(*^+^) = ^ c„ < oo. If 
i',^'k+i) — f{'^k+i)f{'^'k+i) we can bound the l.h.s. of (|8^ by the sum 

j dxfe+idx'^_^i j('')(xfc;x'fc) ((5^i(4+i - Xk+i) - S{x'k+i - Xk+i)) 5p^{xj - Xk+i) f {xk+i)J {x'^^i 
j dxfc+idxfc+i J^^\xk;x'k) {Sp^{xj - Xk+i) - S{xj - Xk+i)) S{x'k+i - Xfe+i)/(xfc+i)7(xfe+i) 



(8.3) 



We next bound the first term by 

dxfc+idx'fc_^i J^''\xk;x^)5/3^{xj - Xk+i){5/3-,ix'k+^ - Xk+i) - S{x'k+^ - x^+i)) f {xk+i)J 
< ll'^^'^^lloo j dxk+idxf,5i3^{xj - Xk+i)\f{xk+i)\ 

X /(x'fc,Xfc+i)- / d4+i5^,(a;fc+i - 4+i)/(x';,,x'fc+i) 



(8.4) 

A slight generalization of a standard Poincare-type inequality (see, e.g. Lemma 7.16 in [7]) yields 
that 



/(xfc,Xfc+i)- / dx'k+iSf3{xk+i - x'k+i) f{xk,x'k+i) <C 



\^k+if{^'k,xk+i + y)\ 
\y\' 



dy (8.5) 
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holds for any x'^ and Xk+i- Inserting this inequahty on the r.h.s. of ()8.4() and applying a Schwarz 
inequality we get 

j dxfc+idx'fc_^i j('')(xfc;x'fc) (5/3i(4+i " ^k+i) - (^(4+1 " Xk+i)) 5p^{xj - Xk+i) f {xk+i)J (x'^.^^] 
< ||j(^)||oo j dxfc+idx',dy 



</3i) 



Xj Xk+l 



+ 117(^^)1100 / dx',drEfc+idy^^M|A)|Vfc+i/(x'„Xfc+i+y)p 

"/ I I 



where in the first term we integrated over dx^ and over dy (the integration over dy gives the factor 
while in the second term we integrated over dx^ (in particular over dxj). To control the first 
term on the r.h.s. of the last equation we apply (|5.4p . In the second term we shift the Xk+i variable 
and then we integrate over dy. We find, using that U^''~^^\xk+i,x'i^^j^) = /(xfc+i)/(x'^_|_]^), 

j dxfc+idx'fc+i j('')(xfc;x'fc) (5/3i(4+i - ^^fc+i) - S{x'k+i - Xfc+i)) Sp^ixj - Xk+i) f {xk+i)J {x'^+i) 

< Cp^'^W^ (3i Tr (1 - A,+i)(l - A,)C/('+') 

(8.6) 

uniformly in /?2. In order to control the second term on the r.h.s. of (|8.3() . we use that 

j dxfc+idx';,+i j('')(xfc;Xfc)(5(xfc+i - Xk+i) (Sp^ixj - Xk+i) - 6{xj - Xk+i)) f (xk+i)! (x,^^^) 

= J dxfcdx'fcdxfc+i J^^\xk;x'f.) i6^.^{xj - x^+i) - 5{xj - Xk+i)) f{xk,Xk+i)f{x'i„Xk+i) 

and then we apply Lemma 18.31 below, keeping all variables fixed, apart from Xj and x^+i- This 
completes the proof of the proposition. □ 

Lemma 8.3. Assume 5p is defined as in Provosition \8.1l Then we have 
J dxdzdx' J(x, x'){6i3{x — z) — 5{x — z))f{x, z)f{x', z) 

< II J|U + ||V, Jlloo ] I J dxdzJix, z){l - A,)(l - A,) fix, 

Proof. From the version (|8.5|) of the Poincare inequality, we can bound the left side of 1)8. 7() by 



(S7) 



\x — z\ 



■\fix\z)\ 



<C dzdx' / d2;l(|x-z| < 2/3) 



\V,{J{x,x')f{x,z)]\' , _l |/(x^z)p 

^ i To ~r ^ i To 



\X — Z\ 



\X — Z\ 



In the first term we drop the restriction l(|x — z| < 2/3) and apply the Hardy- type inequality (|5.2() 
to the z variable on the domain A. In the second term we perform the dx integration. Therefore the 
l.h.s. of (|8.7|) is controlled by 



C J dzdx' J dx K|V^[J(x,x')V^/(x,z)]p + k|V^[J(x,x')/(x,z)]|2 + i^-^f3\f{x' , z) 



(8.9) 
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The first two terms of the last expression are bounded by 

l2 



Ck\\J\\oo + IIV^JI 



f{x,z){l - A,){1 - A,)f{x,z)dxdz 



while the last term is bounded by C/3k ^ll/lli- Optimizing the choice of k, we obtain ()8.7|) . □ 
As a corollary of the proof of this lemma, we can prove that for all < /3 < 1, 

dxdz6p{x-z)\f{x,z)\^ <C j dxdz7(x,z)(l-A,)(l-A,)/(x,z) 

with the constant C independent of (3. This reproves the second part of Lemma 15.21 (originally 
proven in Lemma 5.3 of ^) for the case of radially symmetric potential U and integration over a 
finite volume A. 



9 Proof of the Main Theorem 

The strategy for the proof of Theorem 12.11 is as follows. First, in Section [9.11 we derive a hierarchy 

(k) 

of equation for the time evolution of the densities Uj^^^ for large, but finite A^. 

Using this Hierarchy we prove in Section PO that the sequence U^.t = {t^lfllfcLi is equicontinuous 
in t with respect to the weak* topology of H^^^ . This result will be used to prove that [/tv,* is compact 
in the space C([0, T], H^^^) with respect to the topology induced by the metric p (defined in (|2.1(ij) ). 

In Section f9.3l we prove then that any limit point of U^^t in C{\f^,T], H^^^) (w.r.t. the metric p) 

is also a limit point of the sequence F^v,* = {T/^^I^i (and that any limit point of Tj^^t is also a limit 
point of Ui^^t)- This result implies the compactness of Fjv j. Also part ii) and iii) follow immediately 
from the fact that the limit points of Un^i and Ftv,* coincide, and from the results of Lemma l7. II fin 
particular, the bound (|7.2j) ) and Lemma E31 Having part ii), part iv) of Theorem 12.11 follows easily 
from Proposition 18.11 

Finally, in Section 19. 4[ we complete the proof of Theorem 12. H by proving part v) : here we will 
start from the hierarchy we are going to derive in Section [9.11 and we will take the limit N ^ oo, 
using again Proposition 18.11 

9.1 Convergence to a Regularized Gross-Pitaevskii Hierarchy 

In this section we begin our analysis of the hierarchy of equations governing the time evolution of 

(k) 

the densities Uj^^, for large A^. Here and henceforth we use the pairing 

(J('=),[/W) = I dxfcdx',/'^(xfc;x',)^W(x,;x',). 

The following lemma will be used in Sections 19.21 and 19.41 in order to prove the equicontinuity of 
U^^\ with respect to the weak* topology of H^^ and to prove that any limit point of T]\f^t satisfies 
the Gross-Pitaevskii Hierarchy (|2.241) . 

Proposition 9.1. Suppose J^^'^ G W^''^{h!^ x A'^). For < /? < 1 we choose a radially symmetric 
function 5p G C^(M^) such that < 5p{y) < C(3~'^x{\y\ < and J dy5p{y) = 1. Then, for any 
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t > 0, 



-STTiaoV] / / dxfcdx'fcj('')(xfc;x'fc) 



(9.1) 



X / da;fc+i((^^(xj - Xfc+i) - - Xfc+i))C/^^^''(xfc,a;fc+i;x'fc,a;fc+i) 



as ^ cxD. 

Proof. Put = (f>i\f^t{^)- From 
(see (jnSl) and (jHIZl)) it follows that 



idt(pt = L(j)t + B4>t 



dx^_fc(VFW)=^ 



{LcPtW^ - {L'cP'^)c^t )+{B- B')M't 



(9.2) 



where the superscript ' means that the coordinates xi, .., x^, are replaced by x'^, .., x'^. Similarly, we 
will use := '^x'^- In most cases we will not write out the arguments of all functions fully, but 
it is understood that functions (pt, B, etc. with unspecified arguments depend on the variables 
xi,X2, ■ ■ ■ ,X]\f and their primed versions, (/>^, B', Q' , etc. depend on x'i,X2, . . . , x'f^,Xk+i, ■ ■ ■ ,xn- 
Next we note that for any fixed k 



m=l V / 

k N N / \ 

= - ^ - 2 (V„ log Vm- Yl ^™ + 2Vm(log W^''^)Vm J 

■m=l m=l m=k+l \ / 

since W^^'^ is independent of the first k variables. The contribution of the last term Ylm + 
2Vm(log VF['^l)Vm cancels the analogous contribution from L' in (|9.2j) by self adjointness on the 
space of the last N — k variables. Moreover we use the estimate (see 1)6. 

B-B' = {qjm - q'jj + - + 0(exp(-cr^)) 

(with the summation convention). Thus, from ()9.2p we get 



k 



k ^ r 

:J](-A,+A;.)C/i'l(x,;x',)+ J2 J d^M-k{W^'^?{qjm-q^ 

j=l m,j=l 



Vmlog 



w 



(V 



m(Pt)(Pt 



[ d^N-kiW^"^? 

m=l 



V'^log 



w 



(v: 



(9.3) 
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We can rewrite this equation in integral form 



N t 
.,,=1 -^0 



N 



+ 2i V / ds / dx^_,(W^W)2 

^) (Vm</'.)'A'. - (v:„ log j (v;„</>'j,/.. 



(9.4) 



Only the terms in the first three lines of last equation survive in the limit N ^ oo. The other terms 
vanish as N ^ oo, after they are tested against a function J^^\xk; x'f^) £ W^'°°{A^ x A'^): this is 
proven in Lemma llU. II and Lemma llU.2l in Section [TUl Thus we are left with 



j=i Jo J 



N,s 



N t 

V ds J dxkdx'i,dxN-kJ^''\W^^^)\qjm-q'jm)Ms' + toil) 



(9.5) 



j,m=l 



as N ^ oo. 

Next we note that qjm — Qjm = if j, m > k. On the other hand if j, m < k, then we first use a 
Schwarz inequality to separate (j) and (/>', then we use (|5.3j) with U = xil^m — Xj\ < 2ii) and we get 



j,in<k 



E / dXfcdx',dx^_,|j('=)|(H^W)2g„,| 



^ \<P's\ 



<a£r'll«^^'^lloo I dxfcdx',dx;v-fc(^^['l)'x(km-x,| <24)(|</'s|' + |<^ll') 

j,m<k 

< CkaefejWJ^^^Woo J2 / dxfcdx;dx,v-fe(VFW)2 (jVjcPsl^ + 

j dxTy2(|V^.^,|2 + |^^|2)+ j dx' {w'f\<p'. 



(9.6) 



In the last line we used that the volume is finite and that W^''^ < C^W' (by ()('.1|) ). Since a ^ £i by 
assumption, the r.h.s. of 1)9. 6p vanishes in the limit N — > oo. 

The argument above implies that the factor qjm — q'jm only gives an important contribution to 
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H9.5|) j<k<m or m<k<j. Using the permutation symmetry of the wave function we get 

k ... 

N,s 



-2i(iV-A;)V ds dxfcdx'fcdx7v-fe J^''^(Vr['=l)2 
j=iJo J ^ 

X {q{xj -Xk+i) -q{x'j - Xk+i))(j)s{^k,^N-k)4's{^'k^^N-k) + t o{l) 

for — > CO. 

In the next step we replace T^t^l by To this end we use that 

N 

m=k+2 

N i-1 N 

+ E(Gf^-G^") n Gt'^ n 

i=k+2 m=k+2 m=i+l 



and estimate the effect of this difference in ()9.7p . For the first term in 1)9.8^ we use 



r>fc+l 

and its contribution to (|9.7|) for any fixed j < k can be bounded by 



(iV-A;)||J«|U E /dxfcdx'.dxjv-fcu^fc+vi^l?!,, ( H |<A1 

r>fc+l \m=fc+2 / 

<Ca^r'(^-^)ll«^^'^lloo E [ d^kd^'kd^N-kWk+i,rFl'l, ( J] G'w'j 

r>k+l \m=k+2 / 

xxi\xj-Xk+i\<2ei)(\^s\^ + \^'f) 
<CaqHN -k)\\J^^^\U E I d^kd^'kd^N-kWk+i,rFl%^( J] 



rj[k] 

2 / 



(9.7) 



(9.8) 



where we used Xj.k+i ^ Cl\Xj^k+i and we apphed the usual Hardy inequality (jS-lf) for the variable 
Xj. Next we use w < Cah\ r ^ CFk+i,r and we estimate Ilm Cm < C^^+^VF^ (following from 
Lemma IC.1|) . so that we can use \b.7\ from Lemma 15.31 with respect to the variable Xk+i- We find, 
for any fixed j, 



'Mi\\4>s\ \4>s\ 



(N-k) J2 [ d^kd^N^kWk+i,rFl% ( n G^lkM 

r>k+l \m=k+2 J 

< Caq\N - k) ja^i j dxdx'fcdiy2pVfc+iV,-(/>,|2 + \Vk+i(^s? + iV^+i^'Pl (9-9) 



+ailr^ j dxdx'^ 
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Finally we apply the estimate W < C^W' for the terms with (p' to obtain the energy norm. Using 
Corollary 16.21 and the finiteness of the volume we find (with aii <^ l"^) 

j<kr>k+l \m=k+2 / 

< Ck{a + ahr'"^) ^0 as ^ oo . 

The other contributions coming from the second term of 1)9. 8|) can be bounded analogously. From 
(|9.7|1 . and since 

it follows that 

{J^'\ui;l) = {J^'\ui;l)-^f;^ J ds J d^,d^u'^'\^k;^'k){-^j + ^jWi^li^kM 



k ,t 



- 2i{N - A:) V / ds [ dxfcdx',dx,+i J^'H^k; ^'k) (9-10) 

X {q{xj - Xk+i) - q{x'j - Xk+i))U^^^^\-xik, Xk+i;:^.^, Xk+i) 
+ to(l) . 

Note that the function q depends on N, and that Nq(x) approaches 47rao times a Dirac delta function 
as — > oo fLemma |A.2|) . Using Proposition |H3 (with /3i = 0, that is (5/3^(x) = 5{x), and P2 = 3/2£i) 
we can replace Nq{x) by 47rao times a ^-function. We can also replace Nq(x) by 47rao times a 
smoothed version of the (5-function at some fixed length scale /3. More precisely, we choose a radially 
symmetric function 613 G Co°(IR^) such that < Sf3{x) < C(3~^x{\x\ ^ P) ^^d f dx6i3{x) = 1. By 
the assumption J^'^) E W^'°°{A'' x A'^) and the a-priori bound from Lemma l7. II we have 

j ds j dxfcdx'^,dxfc+i j('')(xfc; x'^) [(A^ - k)q{xj - Xk+i) - AnaoSfsixj - Xk+i)] 

X f/(,'+^)(xfc,x,,+i;x'„Xfc+i)| < Ctiel^^ + P'/^) , 

for some constant C depending on k and J^'^'^ but independent of A^ and /3. Here we used Lemma 
18.31 twice, once for the function q with lengthscale ii fLemma IA.2|) and once for the function 6p. 
From 1)9. 1U() we find 



k t 

{J^'\U^n\) - f /"dx,dxU('=nx.;xl)(-A,+A;.)[/^i(x,;x',) 

k pt p 

-SniaoY] / / dx^dx';, j('')(xfc; x'^) 

X j dxA:+l(5/3(Xj - Xfc+l) - 5/3(x^- - Xfc+l))^7^^^^(xfc,Xfc+i;x';,,XA;+i) 

+ tOi(3^/^) + to{l) 

as A^ — > 00. □ 
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9.2 Compactness of U, 



Using Proposition 19.11 we can now prove the equicontinuity of Uj^ [ in t with respect to the metric p 

on C([0, T], /fi'^^), and thus the compactness of UN,t- 
We recall the bound 



\\Un\\ 



C([0,T] 



sup Y.^-^\\U^^lh<l 

tG[0,T] 



(9.11) 



for some sufficiently large v > 1 (Lemma 17. If) and we recall the definition of the metric p from 
1)2. 16(1 . In order to prove the compactness of the sequence UN,t with respect to the topology in- 
duced on C([0, T], ffi'^^) by the metric p, it is enough, by the Arzela-Ascoli Theorem, to prove the 
equicontinuity of the sequence UN,t- 



Lemma 9.2. The sequence of families of density matrices Un,i 



1,2, 



satisfying 



jy.ll]) is equicontinuous on with respect to the metric p if and only if for every fixed k > 1, for 
arbitrary J^^^ £ ]/l/i.°°(A'^ x A'^) and for every e > there exists a 6 > such that 



< £ 



(9.12) 



whenever \t — s\ < 5. 



Proof. Equicontinuity in the metric p means that, for any e > there exists 5 > (independent of 
A^), such that 



p{UN,t, Un,! 



u 



(fc), 



k>l 



< e 



(9.13) 



whenever |t — s| < 5. Recall that Jj = {Jj }fc>i, for j > 1 was chosen as a dense countable 

subset of the unit ball of H^]^\ Using the uniform bound (|9.11j) . one can approximate any given 

J = {J^^\ J^'^\ . . .) G by an appropriate finite linear combinations of Jj and one can easily 

prove that 1^1^ imphes (|^TT^ . 

On the other hand it is clear, by a standard approximation argument (and because W^''^{h}' x A'^) 
is dense in L^(A'' x A'')), that (I^TT^ for all J^") e W^'°°{A^ x A'') implies the same bound for all 
jik) g j2(^j^k j^ky pj,Q^g ^YiaX (limi) for all j(^) G L'^{A^ x A'=) imphes (Pmi) one can proceed 
as follow. Given e > 0, we note that 



^2-lj:(jf,<, 

j>m k>l 



u 



N,si 



< 2 ^ 2"^' < e/3 

j>m 



if we choose m sufficiently large. Hence 



(k) rr(k) 
' ^ N,t 



U 



N,s/ 



k>l 



<./3+ j;2-^'|^(jf 

j<m 



(fe) j^{k) 
N,t 



u 



N,s/ 



(9.14) 



fc>i 



Moreover we note that 



j:2-|j;(jf),f/S 

jr<m k>p 



u 



N,s/ 



<2 sup ||;7iv(t)|| 

te[o,T] 

< e/3 



H 



.)5]2-^-Lp^'ii'^ri 

" j<m V>P 
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til (fc)ll 

if we choose p large enough, depending on m, because linifc^oo Ik) ||a;,+ = for every j. From 



■ 141) we therefore have 



5: 2-^'| 5:(jf ), Ui;i - < 2e/3 + 2-| ^.iJ^'^U^ " <) 

j>l k>l j<™ fc<p 

Now, for every j < m and k < p, we can find > such that 

if |t — s| < 6jk- Thus, for \t — s\ < 5 := mm{6jk ■ j < ni, k < p}, we have 

5:2-1 ^^(Jf <e. 

j>l k>l 

This proves that (jlTT^ imphes ^TH^ . 



(9.15) 



□ 



Lemma 9.3. The sequence UN,t = {U^\}^^i G C{[0,T], H^^) is equicontinuous in t with respect to 
the metric p (defined in (|2.15|) ). In particular, by the Arzela-Ascoli Theorem, the sequence U^^t is 
compact in C{[0,T], H^^) with respect to p. 

Proof. We prove (|9.12|) for ah J*^'^) G VF^'°°(A'^ x A'^). For such J*^'^) we can apply Proposition 19.11 
and we find 



'N,s/ 



< E / 'd. (J('=),(-A,+A;.)f/, 

+ 2ao ^ r ds / dxfcdx';, J^'^^Cxfc; x'J 
1=1 ■'^i 



(9.16) 



Next we bound 



and use 



+ |ti-t2|0(/3'/')+o(|ti-t2|) . 



(jW,(-A, + A;.)C/if]) < Trj('=^A,[/(,') + TrJ^C/J^A,- 



q.Tj-W q. 



(9.17) 



Tr j(^) Aj- f/J^*;^^ <\\S-^j'^^^Sj\\\\Sj^ Aj SJ^\\ Tr 
Using that J^'^^ G VF1''»(A'^ x A'^) and the finiteness of the volume of A we obtain 

J dxfcdx';, |VjJ(xa;;x'^,)P < "O- 

It follows that S'jj(*^) is a Hilbert-Schmidt operator, and thus compact. This implies in particular 
that ^-VWS'j is a bounded operator, and thus, by Lemma 17.11 the r.h.s. of 1)9. 17|) is bounded. 

Analogously, it also follows that \Tt J^^'^ U^\Aj\ is bounded, uniformly in and s. 
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As for the second term on the r.h.s. of ()9.16|) . we use that for fixed /3, 5p is bounded, J^'^^(xjt; x'^ 



is bounded and Tr U^^^^ < 1: we obtain 



j dxfcdxfcdxfc+i J'^''Hxfc;x'^)((5/3(xj - Xk+i) - Spix'^ - XA;+i))C/^^^^(xfc, x^+i; x'^, x^+i) 



< C 
(9.18) 



for a constant C, independent of N and of s. It follows from ()9.16|) that, for any fixed J^^^ G 

(J«,C/«^)-(J«,C/«^) <C\t,-t,\ 

for a constant C which depends on k and on J^^\ but is independent of N and of ti and t2- This 
implies, by Lemma 19.21 that 11^(1) is compact w.r.t. the topology induced by the metric p on 
C{[0,T],H^_^^). 



□ 



9.3 Compactness of F 



The aim of this section is to prove part i) of Theorem 12.11 stating the compactness of the sequence 
T^^t = {7{f,l}f=i in C([0,r],5i''^) with respect to the metric p. 
First of all, we note that, for any v > 2, 

l|r7V,i||^M=E^"'ll^Sll2<l 
fc>l 

because ||7^^^||2 < 1 for every k. Thus Fiv,t € B^^\ 

In order to prove the compactness of Fat t we use that, by Lemma 19.31 the sequence UN,t is 
compact. It only remains to prove that limit points of UN^t are also limit points of F^v,*. This is the 
aim of the next two lemmas. 

Recall that 

are the marginal densities corresponding to ipN,t{x) = W{x)(j)]\f,t{x), while 

J7^j(xfe;x';,) = j dxN-kW^''^{xN-kf(pN,tixk,^N-k)(pN,ti^'k7^N~k)- 

Recall that IF^ = IF(1-'=) denotes the wave function W after removing its dependence on x^ = 
(xi, ...,Xk) (see (1131)). 

Lemma 9.4. Assume that Fjv^j = {'y^]^\}k=i o.f^d UN,t = {U^\}^^i are defined as above and that the 
assumptions of Theorem A2. 1\ are satisfied. Then we have, for every fixed k > 1 and t G [0, T], 



dxfcdx'fc 



7^j(xfc;x'fc) - U^^{{xk;x[ 



as N ^ oo. 



Proof. In this proof k is considered fixed and all constants may depend on it. We have 

dxfcdXfc 7^](xfc;x'^) - ;7^^j(xfc;Xfc) < j dxkdx',^dxN_k\(t>N,t{^k,^N-k)\\4>N,t{^'k^^N-k 

X W{xk,XN-k)WixlxN-k) - PFW(x^-fc)' 



42 



Using the shorthand notation cj) = (j)iy^t{^ki^N-k), 4>' = (kN,t{^'k-:'^N~k)i and W = VF(xfc, x^v-fc) 
W = VF(x'^,X7v-fc), and VK^ = W^^\^N-k), we need to bound 



\WW' -{W^'''^fUU'\< / \W -W^^\W'\(t)\\4>'\+ / VFW|W^W|(/>||(/>'| 



(9.19) 



Using W = nf=i g]'^ and = nf=fc+i(Gf V^', we get 

k N 



k N 
j=l ' j=k+l 



N / m-l \ N 

m=fe+l \ j=fc+l / \ j=m+l 



Thus, the first term on the r.h.s. of (|9.19j) can be bounded by 



N 



N 

+ E 

m=k+l ' 



I \W -W^^\w'm<l^'\ < I \i-\{g]'^\ n Gf\w'\ct^\W\ 

y \j=k+l ) \j=m+l j 

Applying Schwarz inequahty, with some a > which will be specified later on, we find 

j\W-W^^\w'\^\\ct>'\<aj\l-\{G]'^n n G,j|</.|2 

j=i v=fc+i / 

N „ / m-l \ 1 ( ^ 



r)i=fc+l \i=A:+l 
k 



i j=ni+l 



< 



TV 



m=k+l ' 



(9.20) 



Next we note that, since < G,- < 1, we have 



j=l j<k m 

Using 1)5. 7|) summed up for all 1 < j < A; with the help of (|B.2|) . the first term on the r.h.s. of 1)9. 2U() 
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can be bounded by 



i=l j<k m 

<ahY^ j W'^\Vj(t)\^ + allr^k j W 



for iV — > oo and for some 6 > (because ii <C i^^'^). As for the second term on the r.h.s. of H9.2U() 
we note that, since Gj > ci > (pointwise, for large enough), we have for each fixed m 

^ ^ ^ Wran Pmn ~l" ^ ^ '^mn \ Fmn ^rrm I 



n<k 



n<fc n>kT<.k 

using HC.4|) . Thus, using Lemma IB. II and the fact that tt; < X; we find, 

m>k n\,n2<k 



<n>k m>kn\,n2<k 

+ Cl-^' E E j W^W^mn,WmnAF^X(^mr,e^rM 

m,ni,n2>k ri,r2<k 

^ ^ E E / ^'^Li^mnl^P 



m>k n<k ' 



+ 



™ u u 'J 



in,n>k r<k ' 

Using the estimates (|5.6|) . Lemma 1^31 and 1)5. 3() with U = 9mr, we obtain 



n<k 

+ cr^%e\\oge\)^a^Yl E / w^'(|v™v,(^p + ^iv^^p) 

m>fe r<fc 

+ C7r2^(£|log£|)2a2^ir3Af^ /" wWVr(l)f + \(l)\^) + 0{£ 



r<k 



< C(^a^ + a^iir^ + ar2^(£| log^l)^ + aiil^^-^'l log^l^) + 0(^^-^) 

as — > oo, for some sufficiently small 6 > (here we use (limi) . that £i < ^3/2^ ^ £4 g^^j ^jjj^^ 
e < 1/10). So choosing a = Ni~^ (for some 6 small enough), we get, from (|9.2()|) . 

|p^_-p^W|^'|<^||<^'| < o(l) 
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for N oo. Analogously we can bound the second term on the r.h.s. of H9.19p . □ 

Lemma 9.5. For any increasing subsequence Nj, the subsequence U]\f.^t converges if and only ifT^.^t 

converges (the convergence is in C([0, T], i/i*^^) with respect to the metric p). Moreover the limits 
coincide. 

Proof. Suppose that, for a given subsequence Nj, U^j^t C^oo,t = {f^i^^t}fc>i as j oo, with respect 
to the metric p. Then we prove that '^Nj,t f^oo,t w.r.t. p for j oo. Since '^Nj,t G B^-'^ (the unit 
ball of H^^) it is enough to prove that for every fixed A; > 1 and t G [0, T], and for all J^^^ from any 
dense subset of Lp'{K^ x A'^), 

j dxfcdx', J(')(xfc;x',) (7j^yxfc;x',) - [/£i(xfc;x',)) ^ 

as j ^ oo. Assume now that j'^'^^ G W^'°^{A'^ x A*^) (which is a dense subset of L'^{A'^ x A'^)). Then 
we have 



/ 



dxfcdx';, j('')(xfc;Xfe)(7)J^^^j(xfc;x'^,) - ^/^^^(xfc; x'^,)) 

< j dxfcdxfc 7{f^yxA,.;x'/,,) - [/^^.yxfc;x'fc) 

+ j dxfcdx'fc j(*^)(xfc;x'fc) (^[/Jf^.^Xfe; x'^) - C/^_^i(xfc; x'^; 

The second term converges to zero, as j — > oo, because we assumed that UNj,t Uoo,t w.r.t. the 
metric /? as j — > oo (and because UNj,t S B^^). The first term converges to zero as j ^ oo, by 
Lemma This proves that ^Nj,t — > ^oo,t w.r.t. p as j — > oo. Analogously one can prove that, if 
rAfj,t Too,*, then also UNj,t ^oo,t- □ 

The last lemma, together with Lemma 19.31 implies that the sequence Ti^^t is compact, and 
completes the proof of part i) of Theorem 12. II 

Part ii) of Theorem 12. II follows from Lemma 19.51 and from ()7.2() . 

Part iii) of Theorem 12.11 follows on the other hand by Lemmas 19.51 and 17.31 

Part iv) follows by the remark after Proposition 18.11 using that Foo,* satisfies the bound from 
part ii) of Theorem \2l\ Tr (1 - Ai)(l - Aj)-f^^\ < , for all i / j. 

9.4 Convergence to the Gross-Pitaevskii Hierarchy 

In this section we prove that the limit point Uoo,t satisfies the Gross-Pitaevskii Hierarchy, in the 
sense of fTM . 

Proof of part v) of Theorem, \2.1\ Using the assumption j'^'^-* G VF^'°°(A'^ x A'^) we can apply Propo- 
sition IHUI We find 

k 



(JW, = (jW,[/(fi)-z^ fds j dxfcdx',jW(x,;xl)(-A,+A;.)C/i'^,(x,;xl) 

- Smap j dsy dxfcdx'fcj('')(xfc;x'fc) {^.21) 



+ tOiP^/^)+to{l) 
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as ^ oo. 

By passing to a subsequence and by Lemma 19.51 we can assume here that UN,t — ^oo,t = 
{7oo\}fc>i with respect to the topology induced by the metric p on C([0, T], H^^). This in particular 

(k) (k) 

implies that Uj^^ — > j for every k > 1 and for every t £ [0,T] w.r.t. the weak topology of 
^2(^fc X A'^). Since J^'^) e W^'^{A'' x A'^) and |A| < 00, we have jW(xfc;x'^) e L^{A'' x A'^). This 
means that 

{J^'\ui,'l-jS,)^0 and (j('=),f/i,^l,-7So)-0 

as — > 00. 

Moreover, since jC^) € H^2,oo(^fc ^ ^fc^^ ^^^g^ foUo^g ^hat , j('=)(xfc; x'. ) and A^/ j('')(xfc;x', ) 

are elements of L'^{A!^ x A*"'). This and the fact that Uj^^t ^oo,t w.r.t. the topology induced by the 
metric p, imply that 



sup [ dx,dx',A,j(^')(xfc;x',)([/J,')(xfc;x',)-7g,(x,;x',; 



as — > 00. 

Finally we consider the limit ^ 00 of the last term on the r.h.s. of (|9.2ip . From the proof of 
Proposition 18.11 (see in particular (|8.6j) ). we have 

j dxfcdx'^,dxfc+idx';,+ij('')(xfc;x'fc)5/3(xj - x^+i) 

X {^ri{xk+i — a^fe+i) — ^Xk+i — 2:^+1)) ■'(xfc, Xfc+i; x';,, 

< Cr? II J^IUTr Sk+iSjUl^;^'^SjSk+i 
for some finite constant C independent of s G [0,T], of A^, and of /?. In particular, 

J dxfcdx'fcdxfc+i J^^\xk;xk) {Spixj - Xk+i) - 5p{xj - Xk+i)) [/^^^^^''(x^, Xfc+i; x^, x^+i) 
= j dxfcdx';,dxfc+idxfc+i j('')(xfc; x'/,,) {dpixj - Xk+i) - dpix'j - Xk+i)) 

X 5^{xk+i - Xfc+i)C/Jf+^^(xfc,a;fe+i;Xfc,4+i) + 0(r?) . 

On the r.h.s. of the last equation we can now let A^ — > 00 keeping (3 and rj fixed. By the assumpt 
on J*^^) and by the choice of the functions (5/3 and 5rj, it is easy to see that J('=)(xfc;x'j5/3(a 
Xk+i)^r){xk+i — x'k+i) is an element of L^(A'^+^ x A'^^^) for any fixed (3 and rj. Hence 

j dxfcdx'fcdxfc+idx'^^;^ j('')(xfc; x'^,) {Spixj - xu+i) - Sfsix'j - Xk+i)) S.^{xk+i - x'^+i) 
X (f^i/^^^^(xfc,a:fc+i;x;,Xfc+i) -7(^'+^)(xfc,Xfc+i;x;,x'/,,+i)) ^0 



,10ns 

Xj — 
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for N — > oo, uniformly in s. So, after taking the limit N oo, H9.21() becomes 



X {di3{xj - Xk+i) - 5p{x'j - Xfc+i))7^+^)(xfe 
+ 0(/3V2) + o(r?) 



/ 



4+1) 



for any fixed t and k. Next we apply Proposition 18. II to replace — x'^j^^) by (^(x^+i — 

and 5i3{xj - Xk+i) (respectively, 5p{x'^ -Xk+i)) by 5{xj - Xk+i) (respectively, by 5{x'- -Xk+i)). The 

error here is of order + rj. Hence, letting r/ ^ and /3 ^ we find 

{J'^'\lt\) = {j'^'\l%)-^Y. f f dxfedxlj('=)(xfe;x',)(-A,+A;.)7gs(xfc;x',) 

-STrzaoV] / ds / dxfcdx'^dxfe+i /''^(xfe; x';,) 
,=1^0 7 

X -Xk+i) -S{x'j -Xfc+i))7^+^)(xi,...,Xfc+i;xi,...,Xfc,a;fc+i) . 

rn 



10 Control of Some Error Terms 

In this section we use the notation 



V{cP) = jdxW"^ (|ViV2</>|2 + N-'^lVjcPl'^ + |Vi,^|2 + |0|2) . 

By Corollary O ^(0) < C' if {W(I),HW())) < CiN and {W(I),H^W(I)) < CgiV^, 

Lemma 10.1. Assume a <^ ii <^ £ <^ 1, aii <^ ii <^ Then for any fixed k and any 

j{k) ^ |yO,oo(yYfc X A'^')^ and for every wave function (j) symmetric w.r.t. permutation, we have 

I dXfcdx'fcdx^.fc j('=)(Xfc,x'fc)(^W)2(J^ - f]')</'(Xfc,X^_fc)^(x';„XjV-fc)| < o{l)V{(P) (10.1) 

as N ^ 00 

Proof. We recall that 



Using the definitions of Up, Qpj and T from the beginning of Section El and the estimates from the 



47 



end of Appendix O we find 



i j<k 



i j<k 

k 

i,j>k r,r'=l 
k 

i,j>k r,r'=l 



(10.2) 



Tliis inequality relies on the fact that if both indices i,j>k then there is a cancellation between Fij 
and Flj. 

Using \w\ < CaliX <^ CaiX and \Vw\ < CaX, it is sufficient to control the first and the third 
terms on the r.h.s of (|10.2j) with |(Vi(;)jj| replaced with CaXij. 

Inserting the first term on the r.h.s. of the last equation into (jlU.lj) . taking the absolute value, 
and estimating J^''^ by its sup-norm, we find that this contribution is bounded by 



i j<k' 



< 



< 



Cai-^YYl j dxfcdx',dx^_fc(TyW)2A,,i^f |c/<(x'„ x;v-fc)| \<t^{^k,^N-k)\ 

Car'YY. I dxfcdx'fcdx^_fc(T^W)2A,,4/' (|</.|2 + |</.f ) 

i j<k 

I I J,' 1 2 



i j<k 



+ ahi-^ j dxfcdx'fc dxjv-fe {W^'^fxijFl. 

+ C^kaiir^ (^j dxVF2|0|2 + j dx'(VF 
= o(l)P(<A) 



^2 I i/|2 



because alil ^ <C 1. Here x' = (x'^, xjv_fc). 

The contribution from the third term on the r.h.s. of ()1U.2() can be bounded as follows. Also here 
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we estimate j'^^^ by its sup-norm. 

i,j>k r<k 



i,j>k r<k ' 

< Car'-' EE/ d^kd^',d^N-k{W^'^fxij{Fl^^)'^%r + \Vj(t>'\') 

i,j>k r<k 

+ Cahr'-' E E/ dx,dx',dx^_,(T^W)2x,,(Fi;V/X (I./)!' + 10? ) • 



i,j>k r<.k 

Next we apply Lemma 15.21 We find 



Cr'-' Y.Y.j d^kd^N~k{W^'''^f\,F^^9„WM 

i,j>k r<k 

^ ^ E E log^l' / dx(H-W)2x., (fW)V' (|V.V,</.p + IV. 



i,j>k r<k 
i,j>kr<k 

X / dx(W^W)2x.,(i^i;V/' + |V.c/>|2 + |V.<^P + |<A|^) 

+ CkJ2 {ai-'~'f\logi\' J dx'(VFW)%(FW)V4|v,.<^f 

i,j>k 



< Ck{(.'-'\ log + 1-^1 log = 0(l)P(</<) . 



We used the assumption <C ^^^^ (and 0<e<l/4). □ 

Lemma 10.2. Assume a <C < ^, < N(.\ < f^/'', and iV^^^ < 1. T/ien /or any fixed k, 
any any wave function (j) symmetric w.r.t. permutations, we have 

< o{l)V{(p) (10.3) 



dxfcdx',dx^_,j(^-)(xfc,x',)(TyW)2 (^V^log:^) (V, 

as N ^ oo. 
Proof. We have 

Using that for each fixed j 

Gj = 1- WjnFjn and cf^ = 1 - Wj^Fj^^xl^ > k) for j > A; 
we have, for each fixed j > k and m, 

|G,- - Gf I < Y.W,nF,n + Y.'^jn\Fjn " F,' J ' 
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and 

+ \ {^w)jm\FjmX{fn < k) . 
We also use that Wjn is supported OH I X j X Yi 

I < 2£i, and 

X{\Xj - Xn\ < 2ii) Fjn Gj = xi\Xj - Xn\ < 2ii) Fjn (1 - WjnFjn) + 0{exp{- cr')) . 

Estimating J^^^ by its sup-norm, we find 

^1 / dXfcdx',dx^_fcj('=)(xfc,x;)(TyW)2 fvrrJog^) (V^(/.)( 

<Y.Y. I dXkdx'kdxN--k{W^''^f{\{Vw)jm\Fjm + Wjn\VmFjn\)\Vr 
m j<k 

+ f dxfcdx',dx;v-fc {W^'^ fw,n\Vm{Fjn " i^J'bl|Vm<^| l-Z-'l ^^^'^^ 

m n,j>k 

+ EE / dxfcdx;dx^_,(VP'['l)'|(Vu;),„||(F,^-FjS)||V™</.||0'| 

m j>A: 

+ O(exp(-cr^))P(0) . 

Applying a Schwarz inequality we can bound the first term on the r.h.s of the last equation by (using 
the summation convention for m) 



'm4)\\<P'\ 



I dxfcdx;dx^_fc(wW)2|(V«;),^|F,,„|V, 

<aa^ / dxfcdx'fcdx7v-fc(wW)2Aj>„F,„|V^0p 
3<k-' 

aa-^Yl [ d^kd^'kd^N-kiW^''^? XjmFjm\(l)'\ 



i<k- 

+ 



To bound the first term on the r.h.s. of the last equation we use Lemma 15.31 to bound the second 
one, we estimate Fj^n ^ and, for each fixed m, we integrate over the variable Xj using that ^ 
doesn't depend on it. We get 



■ mv^||</''| 



^ / dxfcdx'fcdx;v-fe(W^t''l)'|(V?/7),^^|F,>„|V, 

<CaaY, [ dxfcdx'fcdxjv-fe(W^['=l)'xjm^;llViV. 

+ Caaeie-^Y1 I dxfedx;dx;v-fc(VFW)2x,>„F,^Vm</'l' 



j<k' 

+ Caa~^iiN f dx'(TyW)2|(/,'|2 



< Ck {aaN + aahr^N + a'^aliN) V{(P) = o{l)V{(P) 
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where we used a = tl'^ and that £i <C ^^/^ Next we consider the second term on the r.h.s. of (|1U.4() . 
We have, applying first Lemma 15.31 and then Eq. (|5.3() from Lemma 15.21 (using the summation 
convention for the indices m, n) 

[ d^kd^'kd^N-kiW^''^fwjn\ymFjn\\Vmm'\ 



3<k 



< CaiiY, / dxfcdx'fcdxjv_fc(W'''')'A,„|V^F,-„| (|V„</.|2 + 



< 



j<k'^ 

< Ckihr' + £lr^f\iogi\^)v{4>) = o{i)v{^) 

for A'" — > oo, because £i ^ £. 

We consider now the second term on the r.h.s. of (|1().4|) . The bound ()('.5|) imphes that (here we 
are summing over n,j > k and r < k) 

J dXfcx',dx,V-fc {W^'^fwjnlVmiFjn " i^j"' ) 1 1 V„0| 1 0' | 

< Cr'~' I dxfcdx'.dx^.fc {W^'^fw.nFjl'Ojr (ajV^^p + a-i|0f ) + 0(^^-^-^) 

< Ckaalil-'-' j d^{W^^^fxjnFji''\VjVm(t>\^ 
+ Caal\r^-' I dx {W^'^fxjnFji'ejrlVn.cpl^ 

+ Ca-^a^iri-^ I dx'dxfe (H^['=l)2xini^,'i'^,r|V,-0f 

Applying Lemma 15.21 we find 

J dXfcdx^.fc {W^''^fwjn\Vm{Fjn " ) 1 1 V„^</'| | | 

<Caaiir^-'k j dx(VF['=l)2|VjV™0|2 



+ Caa^?r^-^^2|iQg^|2^ y dx(H^W)2(|VjV„,(/>P + A^|V„,(/>P) 

+ cka-^a£lr^-'e^\iogef J dx(w^W)2 (|Vj-(/>f + iv|(/)f) +0(^^-1-") 

< Cfcl log£|^(aa£ir^-^Af2 + aa^^r^-^iV^ + a-^ahl'^-'N 

< CkiaNhr^-^ + a-^^?r ^"^)P(0) < Ckil^^N^/^r^~'V{(t>) = o{l)V{(t>) 

51 



for iV ^ oo. Here we first used tliat lif < (as follows from N'^l^ < 1 and a <. li <. I <^ I), 

then we optimized the choice of a and we used that A^^^ ^ ^^/^ (and that e < 1/8). 

Finally we consider the third term on the r.h.s. of (|lU.4jl . Using (|C.4() we find (using the 
summation convention for m) 

^ / dxfcdx',dx^_fc(wW)2|(Vt/;),^||(F,-„-F]^)||V^0||</.'| 



j>k 



<aY,J2j dxfcdx'fcdx7v-fc(W['=l)^A (fJ2)^/^|v, 



j>k r<k ' 

Next we apply a weighted Schwarz inequality. We find 



^ / dxfcdx',dx^_fc(H^W)2|(V«;),-^||(F,-„-Fj^)||V„0||</.'| 

^"EE/ dXfcdx',dx;v-fc(H^t'l)'^mr.Xim(i^jt!)'/'|V 
j>k r<k 

+ / d^kd^'kd^N-k{W^'^fa^,e^r{Fl'J^)y^\4>' 



j>k r<.k ' 

where we used that aAj{^ < and that < fimj (recall the definition of ajm from (|5.5|l ). In 
the first term on the r.h.s. of the last equation we can sum over the index j (using Lemma |B.1|) and 
then we can apply Lemma 15.21 for the integration over the variable Xr- As for the second term on the 
r.h.s. of the last equation we can first integrate over the variable Xr (using that (j)' is independent of 
Xr) and then we can apply 1)5. 9(1 . We find (using the summation convention for the index m), 

[ dxfcdx;dx;v-fc(H^t'l)'|(Vw;),™||(F,-^-FjS)||V^0| 

<Caf\logi\'^Yl j dx(TyW)2|V^V^(/)|2 

r<k 

+ Ca-^ae^ {log N)Y, [ dx' (VFW)2| Vj V„</.f 
j>k-^ 

< Ckilog N){aNf + a~^l^N)V{cp) = o{l)V{^) 
where we chose a = and we used that N'^f' <C 1. This completes the proof of the lemma. □ 



A Properties of the Two Body Problem 

We set r = \x\ for x G M'^ and for a radial function g{x) on we use the notation g{r) for the 
function r ^ g{x)\\x\=T- Let §'{''') ■= drg{r) be the derivative of this radial function g w.r.t. r. 
Consider the (unique) solution to the zero energy problem 

i)co = 0, [) := -A+ (l/2)y , 

on R^, where lo is radially symmetric and satisfies the condition limr^oo uj{r) = 1. It is known that to 
is always non-negative and by Harnack inequality actually uj > since V is regular. By the maximum 
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principle, the norm of u is also bounded. By writing cj(r) = f{r)/r, the scattering length, a, of 
a potential V is defined by 

a := lim r — /(r) . (A.l) 



When V has a compact support, that is, V{x) = for |x| > R, then /(r) = r — a for r > i?. It is 
known that a is bounded by 

1 

a < ■ 

and also a < R since / > 0. 



a < — / Fdx (A.2) 
ovr 



Lemma A.l. Lei V be a smooth, positive, spherical symmetric function such that suppV C {x £ 
: \x\ < R}, for some R> 0. Let g := (87r)~"^(||y||oo + ||^||li) ctnd denote the scattering length of 
V by a. Let ip be the ground state of the Neumann problem 

{-/\ + \V)^ = (A.3) 

on the sphere of radius L, with the boundary condition 

ip{L) = l, dML) = 0. 

Then, if L is sufficiently large, we have 

i) 

E = 3aL-'\l + 0{l/L)) as L ^ oo . (A.4) 

ii) There is a constant < cq < 1 such that for all \x\ < L we have 

C 

CO < ^{x) < 1, 1 - ifix) < —r (A.5) 

\x\ 

where C is a constant, depending only on the potential. 
Hi) For all \x\ < L we also have the following bounds: 

|V^(x)|< and |vV(x)| <^^^^. (A.6) 

Note that, in ii) and Hi), the constant C is independent of L, if L is large enough. 

Proof, i) We first prove an upper bound for the energy E. We write the zero energy state uj{r) as 
u;(r) = f{r)/r and let 

^(r) := sin(A;/(r))/r. (A. 7) 

Note that 

Q ^^^^ ^ kf'{r)rcos{kf{r)) - sin(fc/(r)) 

Therefore, assuming L > R, ip satisfies Neumann boundary conditions at r = L if and only if 

kL = tan{k{L - a)) (A.8) 

where a is the scattering length defined in ()A.1|) . We define k to be the smallest positive real number 
satisfying equation (jA.Sf) . in particular tp > 0. It is easy to check that there are constants Ci and C2 
such that 

^7l-^l<fc2<5fl + ^V (A.9) 



LH L J - - L^\ L 
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Let i) = -d^ + {l/2)V, then i)f = r^uj = 0, i.e. 

-f"{r) + \v(r)f{T) = . (A.IO) 

In particular, / is linear for r > R and by the normalization limr—,ooUj{r) = 1, we have /'(r) = 1 for 
r > R. Moreover, / satisfies uniform bounds /(r) < C(l + r), ||/'||oo5 ||/"||oo < C* by the boundedness 
of uj and / does not vanish. 

With the help of the identity 

- [sin(A;/(r))]" = k'^f'{rf sm{kf{r)) - kf"{r) cos{kf{r)) , (A.ll) 

we compute 

= ei,^ + k\f' - 1) ( ^}^] \\\- kf" {sin kf ){cos kf) + -y(sin kff . 

The last two terms are supported on < R. Using the equation HA.10|) and the boundedness 
properties of /, we see that the O^k"^) terms cancel in the square bracket: 

^f^^ = k^il,^ + r-^x{r < R)0{k^) (A.12) 

and therefore 

(V,f)V) =A:2(^,V) + 0(A;3) , 

where (•, •) denotes the scalar product on -L^(|x| < L). 
Using the estimate 

smkf{r)>Ckf{r)>Ckr (A.13) 

r > 2R, we also get the lower bound 

(^,^) > 47r / (sinA;/(r))2dr >c[ k'^r^dr = Ck'^L^ . 

Jo J2R 

So we can use ip as a. trial function for the upper bound on the energy 

E<f^<,;' + 0(kL-^') = f,(l+0(L-)). 

This proves the upper bound in (|A.4|) . Before proving the lower bound in (|A.4|) . we prove parts ii) 
and iii) of the lemma. 

ii), iii) We now prove ()A.5|) and HA.6|) . We set m{r) := rip{r) and we rewrite the eigenfunction 
equation ()A.3|) as ^ 

f)m = (-a^ + iy)m = Em . (A.14) 

Since V has compact support, we can explicitly write m(r) for r G {R,L]. Let A := \/^- From the 
boundary conditions on ip, we have 

m{r) = X-^ sin(A(r - L)) + L cos(A(r - L)), R<r<L. (A.15) 

From i) we have A < CL^^^^. This allows us to expand m{r) up to O(A^), we find 

m{r) = r-a + 0{l/L) (A.16) 
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uniformly in r, for r E {R,L]. Using that a < R, this proves that (p{r) = m{r)/r satisfies the 
inequahties in (|A.5|) for r G [R,L\, if L is large enough. The properties (|A.6|) on the interval [i?, L] 
can be easily proved using the explicit formula (|A.15() (and using that a < R). 

We now consider the region r < R. From the Harnack inequality, the ratio between the supremum 
and the infimum of in a given ball is bounded. Since (p satisfies ()A.5|) for r < i2, it follows that 
ip{r) is bounded away from zero for r < R. This proves the lower bound of (^(r) in (|A.5|) . 

To prove the bound fix) < 1 for r < R, consider the ball Bj^ of radius R about the origin. On 
the boundary of this domain ip{x) < 1 — 5 with some 6 > 0, uniformly for all sufficiently large L 
from ()A.16|) . Inside this ball, 

-Aip < -A(p + -Vip = Eip< CL~^ 

since if is bounded by Harnack principle, i.e. — A((/3 + CL~^x'^) < . By maximum principle, 
ip{x) < I — 6 + CL^^R^ < 1 for large enough L for all x G Br. 

Finally we have to prove ()A.6|) for r < R. Since if{0) is bounded, we have m(0) = 0. From 
the equation of m inside r < R, |m"(r)| < Cq. By integrating m' from to r we obtain m'{r) = 
m'(0) + 0{Qr). Integrating once more yields 

m{r) = m'{0)r + O(^r^) , 

hence the bound \ip'\ = 0{g) follows. Differentiating the equation HA.3|) . we obtain \m"'\ < Cg for 
r G [0, R] and integrating three times we obtain 

m" (r) = m"(0) + 0{Qr), m!{r) = m'(0) + m" {fd)r + O(^r^), 

m(r) = m (0)r H h 0{Qr^) , 

and the bound on ip" follows. 

Next we prove the lower bound in (|A.4|) . Given any wave function cf) satisfying the Neumann 
boundary condition at |x| = L, we can write it as <^(x) = g{x)ip{x), where ip is given in ()A.7() . 
and g > satisfies Neumann boundary condition at |x| = L as well. From the identity = 
{i)'ip)g — {A.g)7p — 2V gVif), we have 

/ dx^f)(/)= / dx|V5p-0^ + / (ix\g\^il)^'4) 

J\x\<L J\x\<L J\x\<L 

and from HA.12|) 

/ dx^\)<P > k^\\<pf + f dx\Vg\^ip'^ -Ck^ [ dx^^. 

J\x\<L J\^\<L J\x\<R \A 

Using HA.13() and that / does not vanish, we have •(/'(r) > ck with some positive constant c, thus 
/ dx^\)<p> k^W^'^ + ck^ [ dx\Vg\^-Ck^[ dx^-^ . 

J\x\<L ''\x\<L J\x\<R fI 

Using Hardy's inequality (Lemma 15. 1|) 

[ ^dx<e [ dx \Vgf + ^ / dx \g\^ 



\x\<R I^P J\x\<L L,^ J\x\<L 
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and by -(/"(r) > ck the last term can be controlled by 



■l\x\<L c2L3 y|^|<^ 



This implies the lower bound 

/ d2/#0 > k\l + 0(L-3/2))||^||2 > 1^(1 + o(L-i 

in (EH). □ 



Next we apply last lemma to prove some important properties of the function w{y) and q{y), 
defined in Section ITTl 

Lemma A. 2. Assume Va{x) = {aQ/afV{{aQ/a)x) = N'^V{Nx) (because a = a^/N), where F > 
is a smooth, spherical symmetric potential with scattering length ao and with suppF C {x € M'^ : 
l^^l < Ro}- Let Q = (Svr)"^ (ll^lli + ll^lloo)- Suppose that the functions w{x) and q{x) are defined as 
in ifH)) and ifl3)] . with a <C £i < 1. 

i) There is a constant cq > such that 

Co < 1 - w{x) < 1 (A.17) 

for all X G M^. Moreover 

< Ca^M^^, (A.18) 

\x\ 

for some constant C independent of N . 
a) For X G we have the following bounds on the derivatives of w{x): 

X{\x\ < 3h/2) 



\Vw{x)\ < Ca- 



|V w{x)\ < Cga j^p^^g < Cg 

iiij We have suppt; C {x G M'^ : |x| < 3^i/2}. Moreover 

< g(x) < Caq^xi\x\ < 3^/2) 

and 

|Vg(x)| <C^x(|x| <34/2) 

/or a// X G M^. 
iv) The L} norm of q{x) is given by 

/ g(x)dx = 47ra(l + o(l)) , a^O. 
The constant C in i),ii) and Hi) is independent of N, if N is large enough. 



(A.19) 
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Remark A. 3. Using the functions A(x) and ct(x) introduced in i5.5\) . this lemma in particular proves 
all the bounds in 115. b}) . The bound \Vw\'^ < Caa follows because a{\x\^ + a^)^ < C(|xp + a^)~^ , for 
all X E M?. 

Proof. Let and iPk{x) = 1 — Wk{x) be the lowest Neumann eigenvalue and eigenfunction on the 
ball {|x| < k}, that is 

(-A + -Va{x))ilj^{x) = e^x) 

with the condition that ipnix) = 1 if \x\ = k. tpnix) is then extended to be one, for |x| > k. We 
define <?!)k(x) := ipK{{a/ao)x). Then we have 

(-A + ^V{y))My) = {a/aofeMy) 
for |y| < {ao/a)K =: L, and with 4>k{L) = 1. By Lemma FA. II part i), we have 

(a/ao)^e« = 3aoL-3(l + o(l)) 

and thus 

e^ = 3aK~^{l + o{l)). (A.20) 
From Lemma IA. II part ii) we get immediately 

Co < Mx) < 1 and wJx) < Ca^^^4^^ (A.21) 

\x\ 

where the constants cq and C are independent of k and a (they depend only on the properties of the 
unsealed potential V{x)). Moreover from Lemma lA. II part iii) we find 

|Vw^(x)| <Cga ^^','!;' -"f and \V^w^{x)\ < Cga^^f^^j^^ (A.22) 

where C is independent of k and a. 

From (|A.21|) . taking the average over k w.r.t. the probability measure fi, part i) follows trivially. 
As for part ii), from (|A.22() we have 



\Vwix)\ < I \Vw^{x)\l^idK) < / X(|X| < K)/x(dK) 



(A.23) 



because the measure fi is supported on [ii/2, 3ii/2\. The bound for \V'^w{x)\ in (|A.19|) can be proven 
analogously, using HA.22|) . 

In order to prove iii) and iv) recall that q{x) was defined by 

^^^^ _ J e^x{\x\ < K)ip>^{x)fi{dK) _ f e^x(kl < K)^p^{x)n{dK) 
J ■ilJi^{x)fj,{dK) 1 - w{x) 

From ()A.20|) and since the measure fi is supported on [£i/2, 3^i/2], we get q{x) < Ca£^^. The 
gradient of q can be estimated by 

^ Je^K-^5{\x\ - K)fi{dK) _^ Je^xi\x\ < K)\V^w^{x)\fi{dK) 
^ ^ - 1 - w{x) 1 - w{x) 

1 — w[x) 
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The first term on the r.h.s. of the last equation (where we aheady used the condition ip^ix) = 1 if 
|x| = k) is different from zero only if |x| G [£i/2,3ii/2]. Since the delta function forces k = \x\ we 
find, because of (|A.20() and because 1 — w{x) > cq, 

The second term on the r.h.s. of HA.24|) can be controlled using (|A.20() and HA.22() . The third term 
on the r.h.s. of (|X24l) can be estimated using the bounds ()A.23p and q{x) < Cal^ ^ . This completes 
the proof of iii). 

Finally, to prove the estimate for the norm of we note that 

dx = [dx ^^-^('"'-."^^-j"^^^^"^ = [dx (A.25) 



/ V'K(x)/i(dK) J f ip^{x)fi{dK) 

for a — > 0, because il^nix) = 1 — w^^ix) and 



dx / et,x{\x\ < K)wK{x)fj.{dK) 



<Ca I I dx^^^4^ ) ^(d'^) 



Analogously we can estimate the contribution of Wk{x) in the denominator on the r.h.s. of HA.25|) . 
We get 



dx q{x) = I dx J e^xil^l < K)fi{dK) + a o(l) = / ( / dx xi\x\ < k.) ] n{dK) + a o(l) 
y eKK^/u(dK) + ao(l) 



in 

T 



and thus, with ()A.20j) . 

J dx q{x) = 47ra(l + o(l)) 
for a — > 0, which proves iv). □ 

B Properties of the Triple Cutoff Function Fij 

In this appendix we collect some properties of the function Fij, defined in Section f2. 21 
Lemma B.l (No overlap). For any exponent q > and any fixed i and j ^ f 

\\x^m'F!J\\oo<e-''il'\ (B.l) 

Moreover, for any fixed i, 

Y.^^Fl^^c,. (B.2) 

Proof. Clearly |F| < 1. Suppose there is an overlap between two functions Xij, Xij'^ i-e- for some 
3 / j' we have Xij / and Xij' / 0. Then |xj — Xj|, \xi — Xji\ < £, so 

so any overlap forces Fij to be exponentially small. This proves (jB.!!) and also (|B.2|) . □ 
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We have analogous controls on the derivatives of Fij as well. We define the sum norm of a vector 
X = (xi, . . . xn) G as 



N 



iixii ■.=y: 



k=l 



where |x| is the Euclidean length in M^. Similarly, if A = [Ai, . . . ,A]\j), where Aj is an n-tensor on 
for each j, then 

N 



|A| ■■=J2\^^ 



k=l 

where \Ak\ = \Ak\(^-^3-j0n denotes the tensor norm on (M^)®" derived from the Euclidean norm. In 
particular, if V = (Vi, . . . , Vat) denotes the 3A^ dimensional derivative, and a G N, we have 

fcl k2 ka 



where 



with 



|Vfci . . . Vfc„/|(]R3)®c = sup I (Vfci . . . Vfc^/)(ui (^V2<» ...<^Va) 



j=l iPi 



/3i,.../3c=l i=l 



Vj £ R-^ , \vj\ = 1 



d 



with Vj = {Vji,Vj2,Vj3) 



Lemma B.2 (Control of the derivatives of Fij). Let i <^ 1, i.e. £ < N for some k > 0. For 
sufficiently large N the following pointwise bounds hold: 



i) Let a G N and g > 0. Then for any fixed i,j 



a 7-19/2 



II v"i^ II < c^i-^f: 



(B.3) 



ii) Let K > and recall the definition of O^j from h2.b]) . For any a G N and for k ^ and 
arbitrary m 



Hi ) For any fixed a € N and index k 



(B.4) 
(B.5) 

(B.6) 



iv) For any fixed a, /? G N and index k 



-0-/3 



(B.7) 
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Proof. For i) note that 



for a € N. Then 



using that ze~^ < ce~^/^ with z = cqNl"^ . For higher derivatives the proof is similar: 



1/75/2 



Ff, < c„ 



79/2 



using that 2;"e ^ < with 2; = Mi ^ . 

For the proof of ii) we note that for k ^ i,j the derivative VkFij is smaller than £^~'^ unless 
smaller than K£\ log^|, therefore 



'ki + "kj 
Oki + Gkj 



{hk, + hki) F^. + e""-'-' 
k / i,j . 



Ff + F-'-^ 



2nK~a~e\ 



The proof for higher derivatives is similar. 

As for iii), we use the previous estimate when i,j^k 

i,j^k i,j^k 

We consider the set Sk = {j ■ \xj — Xk\ < -fC£|log£|}. If G Sk, and \xj — Xj'\ < |£|log£|, 
then Fij < exp(— Therefore, apart from exponentially small error, the cardinality of is 
{2K/e)'^, since if there are more j's in the set S^, then \xj — Xj'\ < (e/2)£| log ^|, at least for two 
indices j,j' £ S^- The same argument holds for the i indices, showing that the i summation is only 
over a finite set. This implies that 



i,j^k j&Skd^k if^Ski^k 



Choosing K sufficiently large we obtain the bound HB.6|1 for i,j / k. For the terms with k = i ov 
k = j we use 

MF',^\ < |V°F,«.| < cr-F,f (B.8) 

and the j summation is finite by HB.2|) . 

Finally, the proof of iv) is similar. If k = i, then we have 

jm j 

using (|B.3|) and ()B.2|) and the same estimate holds if A; = j. If i, j 7^ /c, then again the cardinality of 
the index set for i and j is bounded, apart from an exponentially small error. Therefore 



3pK-a-l3-2e 



using (jB 



□ 
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C Local Structure after Particle Removal 



Lemma C.l. There exists a constant Co > such that for any k 

C,'<-^<Co (C.l) 

pointwise for sufficiently large N. Moreover, for any a G N and sufficiently large N we have the 
pointwise estimate 

r\l{kx...ka) 

C,^ < ^ < Co" (C.2) 
uniformly for any family of indices ki, . . . ,ka- 

Proof. Similarly to the proof of (|4.2() . we see that for any fixed k and sufficiently big N > N{k), 

(k) 

all G- are separated away from zero uniformly in i. Therefore it is sufficient to show that 

iy^k i^k my^i,k 

is uniformly bounded for any fixed k. The boundedness of the first term follows from (|B.2|) . As for 
the second term, note that 



1 - e 



i':^<r'[hu + hrnk]Fi^ 



SO 



E E ^^m\F?2-Fim\<2r' ^ X^^h^kF^S 

since Wim < Xim by the support of w. If jxj — > Kl\ \ogl\ for some i, then hik < (-^ and and if 

K IS a, large constant, then £^~^N'^ 0, so this term is negligible even after the summation. Now 
we look at the set 

Sk ■■= {i ■■ i / k, \xk - Xi\ < Ki\ logi\} . 

If i,i' e 5fc, i i' and |xj — < |£|log£| then F^^ < exp(— < exp(— that is expo- 
nentially small. Therefore, modulo exponentially small errors, \Sk\ < {K/{e/2))^, which guarantees 



that the summation over i in (|(L3|) is finite. For each fixed i the summation over m is finite by ()B.2|) . 
The second bound (|C.2|) follows easily by induction. □ 

The same proof immediately gives the following two bounds that are used in the proof. Here k 
is fixed and the constants may depend on k. 



r<k 

and 



E i^i - ^ E(^- + + o(^^-"-^) . (C.5) 

j r<k 

We will often multiply these inequalities by vuim, then on the support of Wim we can use that 9ir and 
6mr are comparable. In particular, we also obtain for each m 

|GW - Gm\ < E E ^rnjF!^]emr + 0(£^-^) . (C.6) 

j>k r<.k 
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D Upper bound for 

The aim of tliis appendix is to prove tliat the assumptions of Theorem 12.11 about the energy distri- 
bution of the initial data tpNfi = WcpNfl are satisfied for a large class of (j)N,o- In particular, in the 
next lemma we prove that {W(j)N,o, H'^Wc/yN^) < CN'^ is satisfied, if the function (pN^ is sufficiently 
smooth. The proof of the inequality {W(j)N,o, HWcpjyfi) < CN for sufficiently smooth (jiNfl is similar 
(but much easier) and therefore omitted. 

Lemma D.l. Assume a ^ li <^ i <^1, aii <^ and l\ <^ ai^^'^ . Suppose moreover that (j) satisfies 

^ 1^(1- Ai){l-Aj){l -Ak){l-Am)(l)< DN^ (D.l) 

i,j,k,Tn 

for some D > 0. Then there is a constant C > such that 

Proof. Following the steps from (|6.15() to (|6.18() (but this time without neglecting the positive con- 
tributions), we find (using the summation convention) 

{W^,H'^W(j)) = j \HW(I)\'^ 

< J W^\ViV,^\^ + C J W^{\ViV,Ge\ + \ViGe\\V,G,\) |Vi(/)||Vj</.| 
+ 2 J W^\B\\Vm^\^ + 2 J W^\VmBM\Vm(l^\ + J W^B^\cl)\^, 



r2 r?2 \ ±\2 



with B = qj^j -\- = Qf^j -\- ^ O ye ^ j . Using Lemmas 16.31 - 1^?^ we get 

j \HW(l)f <C jw^ (|ViV,0|2 + iV|V,0p + iV2|0|2) +C j VF^g. .|v, 
+ C j W'^\ViGi\\VjGi\\Vi(i)\\Vj(^\+ I W^B 
Since W < 1, and by (|D.1() . we find 

J iHWcpl"^ <CN^ + G j \ViGe\\VjGi\\Vi(j)\'^ + G j qij\V 



|2 



+ G j qijqkm\<l)\'^ + C j n 



(D.2) 



2| j,|2 



We start by considering the second term on the r.h.s. of the last equation. Using Lemma iB.ll we 
find 

/|V,G,||V,G,||V.0|^ < C/ (|(V.)„,|^F,„. + |(V„.)„„k.„F„„|V,F.J 

+ Oie-''~') . 



62 



Using \Vw\^ < CA, \Vw\ < CaX, Lemma lR2l Lemma O (with W = 1) and (|^ we find 

\ViGe\\VjGe\\Vicl)f <C J \miFll-\V,<l}\'' + Car^j A™,F„,-|V,0|2 



+ 0{e~'' ), 

which is bounded by CN'^, because a£i£~'^ <C 1, and because of ()D.1|1 . 
Next we consider the third term in (|D.2|) : we obtain 

J %|V™0|2 < Calf J Xij|V^0p 

<CaJ {\ViVjVm^f + N\ViVm(l)\^+N^\Vm^f) 

< CN^ 

using (|5.4|) and (|D.1|) . As for the fourth term on the r.h.s. of i)D.2|) . we have 

<Ca' Yl y"0(l-A.)(l-Ai)(l- Afc)(l- A„ 
+ Ca^l^^ J ^(1 - Ai)(l - Aj)^ 



because of a ^ i'l^'^ and of HD.1|) . 

Finahy we consider the last term on the r.h.s. of ()D.2|) . We use that, by (|().12|1 . < CQkj^im + 
CT^. By Lemma IB . 21 and by the estimate w < alX, \Vw\ < a A, we obtain 



^kj^im\4'\'^ < C j iwsjWrm\AkFsj\\AiFr, 



\{Vw)kj\\{^w)im\\^kFkj\\ViFim\ l^l"' 



<Ca'r' I XkjX^mFl/J'F|i'\cP\' 
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(D.3) 



We note that, if, for example k = i, then j is forced to be equal to m, up to exponentially small 
errors, due to the strong non-overlapping properties of the functions F^j. Hence, up to errors which 
are exponentially small in and using Lemma 15.31 we can estimate the second term on the r.h.s of 
(EH) by 



1/2 ^1/2,. 2 
im It' I 



kj {k,j)^{i,in) 
{k,j)^{i,m) 

< Car\\ogN) j 0(1 - Afc)(l - A,)0 



(k,j)=/={i,m) 



+ Ca^r^ilN^ [ |0|2 , 



which is o{N'^) because a£i <^ i'^. We still have to control the contribution from . Using Lemma 
IB. 21 we have 

y r2|0|2 <C ^ ^|(V^i;),fc||(Vu;) 

mr\WjpWnsFikFjYir \ ^ kFjp \ \ V7- 

Fns\ 

\V„Fik\\VaF^r\\VdF,J\VdFns\ |0r 

J (D.4) 

mr\FikFfYir ~\~ £ WikWrnrFikFiYir) |0| 



which can be bounded by CN"^ in the same way we bounded the r.h.s. of ()D.3|) . This completes the 
proof of the lemma. □ 

Acknowledgements: We would like to thank J. Yngvason for communicating this problem to 
us and for helpful discussions. 
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